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CHAPTER  1 


INTRODUCTION 

1 . 1  BACKGROUND 

Most  river  and  watershed  Improvement  programs  require  knowledge  of 
the  total  sediment  load  likely  to  exist  during  various  flow  conditions. 
To  obtain  this  knowledge  with  sufficient  accuracy,  field  measurements 
are  made  which  include  discharge  and  suspended-sediment  measurements. 
Discharge  measurements  consist  of  depth  and  mean  velocity  measurements 
taken  at  selected  verticals  of  a  channel  cross  section  from  which  the 
total  discharge  can  be  computed.  Suspended-sediment  measurements  are 
usually  made  by  sampling  the  concentration  of  suspended  bed  material  and 
wash  load  with  a  depth-integrated  sampler  at  some  or  all  of  the  verti¬ 
cals  where  mean  velocity  is  measured.  With  stream  discharge  and  sedi¬ 
ment  data,  estimates  of  the  total  sediment  discharge  are  made  using  com¬ 
putational  procedures  such  as  the  modified  Einstein  procedure  (Colby  and 
Hembree  1955),  the  Colby  method  (Colby  1957),  and  the  Toffaletl  method 
(Toffaleti  1969). 

The  acquisition  of  discharge  and  sediment  data  can  require  exten¬ 
sive  field  data  collection  time  and  manpower.  To  reduce  time  and  man¬ 
power  requirements,  automated  sampling  stations  may  be  installed  which 
periodically  record  the  river's  stage  and  collect  point-source 
suspended-sediment  samples.  The  current  practice  for  analyzing  auto¬ 
matically  collected  point-source  sediment  data  Is  to  determine  the 
suspended-sediment  load  using  a  correlation  previously  developed  from 


1 


2 


concurrently  collected  depth-integrated  and  point-aource  suspended- 
aedlaent  samples. 

Lane  and  Kallnske  (1941)  developed  a  computational  procedure  for 
determining  suspended-sediment  load  using  a  point  suspended-sediment 
sample.  Slnc-i  their  research  and  the  development  of  depth- Integrated 
sediment  load  measurement  procedures,  the  determination  of  sediment  load 
from  point-source  suspended-sediment  data  using  theoretically  based  con¬ 
cepts  has  received  little  attention  other  than  by  Brooks  (1963).  The 
procedure  presented  In  this  paper  advances  the  research  efforts  of  Lane, 
Kallnake,  and  Brooks;  eliminates  the  need  for  a  correlation  between 
point-source  and  depth- Integrated  suepended-sedlment  samples;  and  per¬ 
mits  the  direct  computation  of  "total”  sediment  load  (suspended  load 
plus  bed  load)  using  point-source  suspended-sediment  data. 

I . 2  RESEARCH  OBJECTIVES 

The  goal  of  this  research  Is  to  develop  a  theoretically  based  com¬ 
putational  procedure  for  making  total  sediment  load  measurements  using  a 
point-source  suspended-sediment  sample.  Advantages  of  the  proposed  pro¬ 
cedure  over  existing  procedures  Include  the  direct  use  of  automated  sam¬ 
pling  station  data  and  the  reduced  cost  of  data  collection.  To  achieve 
this  goal,  the  following  objectives  have  been  established: 

1.  Develop  e  procedure  for  computing  total  sediment  load  that  uses 
point-source  suspended-sediment  measurements. 

2.  Compare  the  proposed  procedure  to  Colby  end  Hembree’s  depth- 
integrated  sediment  load  measurement  procedure  (Colby  and  Hembree  1955) 
using  published  data. 

3.  Apply  the  procedure  to  newly  acquired  data. 


CHAPTER  2 


LITERATURE  REVIEW 

2.1  NUMERICAL  FORMULATION  OF 

TOTAL  SEDIMENT  LOAD 

Numerical  formulations  that  determine  total  sediment  load  for  a  set 
of  field  discharge  and  sediment  measurements  generally  subdivide  the 
total  load  Into  bed  load  and  suspended  load.  Bed  load  Is  sediment  that 
rolls  or  slides  along  the  channel  bed,  and  suspended  load  Is  sediment 
that  Is  supported  by  the  upward  turbulent  currents  and  remains  suspended 
for  considerable  lengths  of  time.  The  division  betveen  these  two  zones 
of  transport  Is  not  sharply  defined,  but  has  commonly  been  referenced 
and  applied  by  setting  a  bed  load  layer  thickness  of  two  grain  diameters 
as  proposed  by  Einstein  (1950). 

The  total  sediment  discharge  per  unit  channel  width  qT  with  sedi¬ 
ment  load  divided  as  bed  load  and  suspended  load.  Is  determined  by 

d 


where 

a  •  height  above  the  channel  bed  where  suspended  load  begins 
d  •  flow  depth 

V  *  time-averaged  flow  velocity  at  a  distance  y  above  the 
^  channel  bed 


I  V  C  dy  +  q_ 

J  7  y  7  b 
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C  »  tine-averaged  sediment  concentration  at  a  distance  y  above 
^  the  channel  bed 

<lg  ■  bed  load  transport  per  unit  width 
Therefore,  components  which  are  necessary  for  determining  total  sediment 
load  are  the  vertical  velocity  distribution,  the  sediment  concentration 
distribution,  and  the  bed  load  discharge. 

2.1.1  Velocity  Distribution 

The  basis  for  the  velocity  distribution  used  in  most  of  the  known 
sediment  load  measurement  procedures  is  the  Prandtl-von  Karman  velocity 
distribution  equation  (Prandtl  1926,  von  Karman  1930) 


V 

_Z 


2.2 


where 


-  the  shear  velocity,  "VgRS  ,  where  g  is  the  acceleration  of 
gravity,  R  is  the  hydraulic  radius,  and  S  is  the  energy 
grade  line  slope 

k  -  the  universal  von  Karman  constant,  0.4  for  clear  water 
y^  ■  an  integration  constant 

Variations  in  the  application  of  this  relation  stem  from  derivations  of 

the  integration  constant  y^  .  One  of  these  formulations  considers  the 

maximum  flow  velocity  V  to  occur  at  the  channel's  surface  (y  -  d) . 

infix 

Substituting  ^mBX  for  and  d  for  y  in  Equation  2.2  leads  to 


V 

max 


2.3 


A  new  relation  can  now  be  derived  which  does  not  include  y^  by  combin¬ 
ing  Equations  2.2  and  2.3: 
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V  - 
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2.4 


which,  when  Integrated  over  depth  to  determine  the  average  velocity  in 
the  vertical  ,  can  be  expressed  as 


V 

y 


1 

d 


2.5 


Equations  2.4  and  2.5  can  be  combined  and  described  by 


2.6 


2.7 


which  have  both  been  applied  in  sediment  measurement  procedures  (Lane 
and  Kallnake  1941  and  Brooks  1963) . 

Another  variation  In  the  derivation  of  y^  which  has  received  wide 
application  cornea  from  the  research  efforts  of  Keulegan  (1938). 

Keulegan  used  dimensional  reasoning  to  determine  y^  for  hydraulically 
smooth  and  rough  boundaries.  For  a  smooth  wall  he  considered  y^  to 
depend  solely  on  UA  and  v  ,  the  kinematic  viscosity,  such  that 
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m 


2.8 


where  m  is  a  constant.  For  a  hydraulically  rough  wall,  he  established 
a  relation  for  y^  In  the  form 


2.9 


where  ks  Is  the  height  of  the  roughness  elements.  Substituting  either 
Equation  2.8  or  2.9  into  Equation  2.2  leads  to 


V 

JL  . 


1  yD* 

+  -  In  - 

K  V 


2.10 


or 


V 

JL 


.  2.30  .  yD* 
c  +  —  log  — 


2.11 


where  c  Is  a  constant.  Keulegan  applied  the  data  collected  by 
Nlkuradse  (1932,  1933)  to  Equation  2.11  and  formulated  relations  for 
hydraulically  smooth  and  rough  pipe  flow  which  may  be  applied  In  open 
channels.  Keulegan' s  relations  are 


\  y°* 

^  -  5.5  +  5.75  log  ~  2.12 


for  hydraulically  smooth  flow,  and 
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^  -  8.5  +  5.75  log 
*  s 


2.13 


for  hydraulically  rough  flow.  (Keulegan  used  k  »  0.40  when  he  devel¬ 
oped  Equations  2.12  and  2.13.) 

Keulegan  similarly  formulated  relations  for  the  mean  flow  velocity , 
V  In  trapezoidal  cross  sections.  He  determined 


V  RU* 

Jp  -  3.5  +  5.75  log  — 


2.14 


for  hydraulically  smooth  flow  In  pipes  of  circular  cross  section,  and 


V  RU* 

-  3.0  +  5.75  log  — 


2.15 


for  flow  between  smooth  parallel  walls  of  Infinite  extent.  The  mean  of 
these  two  relations  may  be  applied  In  hydraulically  smooth  open  channel 
flow.  The  mean  of  Equations  2.14  and  2.15  Is 


V  RU* 

-  3.25  +  5.75  log  — 


2.16 


For  hydraulically  rough  channels,  Keulegan  formulated 
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J-  -  6.25  +  5.75  log  2.17 

*  s 

using  the  experimental  results  of  Bazin  (1865). 

The  sediment  measurement  procedures  whi :h  use  Reulegan's  velocity 
formulations  apply  them  as  they  were  presented  by  Einstein  (1950). 
Einstein  modified  the  velocity  distribution  Equations  2.12  and  2.13  and 
the  mean  cross-sectional  velocity  Equations  2.16  and  2.17  using 
Nlkuradse's  (1933)  data.  These  modified  relations  are 


V 

JL 


5.75  log 


2.18 


and 


5.75  log 


2.19 


where  x  Is  a  dimensionless  parameter  for  transitioning  between  smooth 

and  rough  flow  and  a  function  of  kg/ 6  which  Is  presented  in  Figure  4 

of  Einstein  (1950).  (6  is  the  thickness  of  the  laminar  sublayer  of  a 

smooth  wall,  11.6v/Uft.) 

2.1.2  Sediment  Concentration 
Distribution _ 

The  differential  equation  for  suspended  sediment  of  uniform  size, 
shape,  and  density  In  two-dimensional,  uniform,  turbulent  flow  is 
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dC 

C  w  +  e  -r^-  -  0  2.20 

y  s  dy 

where  w  is  the  particle  fall  velocity,  eg  is  a  diffusion  coeffi¬ 
cient,  and  the  remaining  variables  are  as  previously  defined.  If  eg 
were  constant  from  the  flow  surface  to  the  channel  bed.  Equation  2.20 
could  be  integrated  to  yield 


C 

JL 
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a 


exp 


-  =-(y  -  a) 
s 


2.21 


where  C  is  a  referenced  suspended-sediment  concentration  at  the 
fl 

distance  a  above  the  streambed.  However,  the  diffusion  coefficient  is 
not  constant  over  depth  and  can  be  shown  to  be 


t  -  6tcU.  £(d  -  y)  2.22 

s  *  a 

where  B  is  a  numerical  constant  relevant  to  diffusion  in  turbulent 
flows  (Vanoni  1975). 

Lane  and  Kalinske  (1941)  assumed  eg  to  be  constant  over  depth  to 
obtain  an  expression  which  was  more  convenient  to  use  than  Equa¬ 
tion  2.22.  They  assumed  the  average  eg  in  any  vertical  section  could 
be  described  as 


2.23 


II 


When  this  expression  is  placed  into  Equation  2.21,  the  sediaent  distri¬ 
bution  equation  becoaes 


C 

_X 

C 


a 


exp 


15w  (y  -  a) 
d 


2.24 


An  expression  for  sediaent  distribution  which  considers  the  diffu¬ 
sion  coefficient  e  to  vary  with  flow  depth  can  be  derived  bv  sub- 

s 

stltuting  Equation  2.22  into  Equation  2.20.  The  derived  expression  is 


c  w  +  J  <d  - 


dC 

dy 


2.25 


Separating  the  variables  and  integrating  Equation  2.25  yields 


2.26 


where 


z 


2.27 


Equation  2.26  was  developed  by  Rouse  (1937)  and  is  used  by  several 
sediment  measurement  procedures  (Colby  and  Hembree  1955,  Brooks  1963, 
Allen  and  Barnes  1975,  Burkhaa  and  Dawdy  1980,  and  Shen  and  Hung  1983). 
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2.1.3  Bed  Load  Transport 

Sediment  In  transport  within  the  bed  layer  is  called  the  bed  load. 
Einstein  (1950)  developed  an  equation  to  compute  this  mode  of  sediment 
transport  which  Is  based  on  the  probability  that  a  particle  will  move 
from  Its  position  on  the  bed  and  another  particle  will  settle  in  Its 
place.  Einstein's  bed  load  equation  is 


-  ♦.v8’3/2”3/2<s.  - 


2.28 


where 

l^q^  •  rate  at  which  a  particle  of  diameter  D  moves  through  a 
unit  width  of  the  bed  layer  per  unit  of  time 

«  the  transport  intensity  for  an  Individual  grain  size  D 
i^  -  the  bed  material  fraction  for  a  given  size  range 
ig  »  the  fraction  of  bed  load  for  a  given  size  range 
pg  »  the  density  of  the  solids 
$g  -  the  specific  gravity  of  the  sediment 
To  apply  Equation  2.28  for  a  given  grain  diameter,  Einstein  assumed  that 
all  the  parameters  on  the  right  side  of  the  equation  were  constant 
except  for  •  The  intensity  of  transport  was  established  to  be 

a  function  of  the  intensity  of  shear  for  individual  grain  diam¬ 
eters.  The  function  which  relates  and  is  presented  by 

Einstein  (1950,  Figure  10)  where  is  determined  by 


pf 


2.29 
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where  £  Is  a  hiding  coefficient  and  a  function  of  D/X  as  given  by 
Einstein  (1950,  Figure  7),  Y  is  a  pressure  correction  parameter  and  a 
function  of  Dgj/6  as  given  by  Einstein  (1950,  Figure  8),  and 
81  -  log  10.6 

8  ■  log  (10.6X/A),  where  X  -  0.77A  if  A/6  >  1.80  and  X 

x  *  1.396  if  A/6  <  1.80  and  A  ■  D,_/x  where  D,_  is  the 

65  65 

grain  size  for  which  65  percent  by  weight  of  the  bed  material 
is  finer 

■  the  density  of  the  flow  medium 
R'  -  the  hydraulic  radius  with  respect  to  grain  roughness 
Equations  2.28  and  2.29  have  experienced  two  significant  modifica¬ 
tions  as  they  have  been  applied  in  total  sediment  measurement  proce¬ 
dures.  The  first  of  these  modifications  was  made  by  Colby  and  Hembree 
(1955)  in  the  development  of  their  procedure  which  has  become  known  as 
the  modified  Einstein  procedure.  Colby  and  Hembree  modified  Equa¬ 
tions  2.28  and  2.29  as  follows:  The  quantity  for  an  individual 

grain  size  D  is  computed  using  the  larger  value  of  the  two  relations 


*  -  l  65  -2* 
*  1  RS 


where  is  the  grain  size  for  which  35  percent  by  weight  of  the  bed 
material  is  finer.  Colby  and  Hembree  arbitrarily  divided  of  Equa¬ 
tion  2.28  by  2  to  fit  their  observed  total  load  data  more  closely. 
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The  second  significant  modification  to  Einstein's  Equations  2.28 
and  2.29  was  made  by  Burkham  and  Dawdy  (1980).  In  Burkham  and  Dawdy 's 
study,  the  relation  for  <|i#  was  redeveloped  and  the  need  to  divide 

of  Equation  2.28  by  2  was  eliminated.  The  basis  for  the  new  rela¬ 

tion  was  Einstein's  (1950)  equation  for  a  uniform  particle  size 


ps  - 


D 

tfs 


2.32 


which  Burkham  and  Dawdy  changed  to 


k 

*.  ■  «•«  S§  2-33 

3 

where  k  replaces  the  uniform  particle  size;  p  is  5.17  slugs/ft  ; 

8  8 

3 

pf  is  1.95  slugs/ft  ;  and  ft  replaces  R'  .  Burkham  and  Dawdy  devel¬ 
oped  an  equation  for  a  nonuniform  bed  material  size  gradation  with  the 
assumption  that  Einstein's  hiding  coefficient  and  lift  correction  param¬ 
eter  can  be  represented  by  the  single  parameter  (C?)^  »  The  resulting 
relation  for  for  a  nonuniform  size  gradation  is 

B  2'34 


where 


8^  •  log  10.6 


2.35 


*  log  10. 6x 


6  +  C  D 
s 


when  D  <  5  +  C  D 
8 


when  D  >  j  +  C  D 
-  s 


Cs  -  -0.62  +  3.12(Dft5  D35)j 


Burkham  and  Dawdy  developed  Equation  2.38  using  data  sets  from  the 

Niobrara  and  the  Middle  Loup  Rivers  of  Nebraska. 

2.2  SEDIMENT  LOAD  MEASUREMENT 
PROCEDURES 

Sediment  discharge  measurement  of  streams  is  necessary  to  determine 
sediment  loads  and  to  develop  or  check  sediment  discharge  relationships. 
If  sediment  were  uniformly  distributed  within  a  channel  cross  section, 
the  measurement  and  calculation  of  sediment  discharge  could  easily  be 
computed  as  the  product  of  a  known  sediment  concentration  and  stream 
discharge.  However,  sediment  Is  not  uniformly  distributed,  and  a  mea¬ 
surement  and  computational  procedure  Is  therefore  required  to  determine 
total  sediment  loads. 

2.2.1  Sediment  Load  Calculations 
using  Point  Sampled  Data 


The  earliest  known  record  of  suspended-sediment  sampling  is  the 


work  of  Gorsse  and  Subuors  In  1808  for  the  Rhone  River  at  Arles,  France 


(US  Inter-Agency  Committee  on  Water  Resources,  Subcommittee  on  Sedi¬ 
mentation  1940).  From  1808  thru  the  1930's,  suspended-sediment  samples 
were  primarily  collected  from  one  or  more  points  in  the  vertical.  Two 
theoretically  based  procedures  which  can  apply  such  samples  for  com¬ 
puting  suspended-sediment  load  have  been  developed  by  Lane  and  Kalinske 
(1941)  and  Brooks  (1963). 

2. 2. 1.1  Lane  and  Kalinske.  Following  the  development  of  turbu¬ 
lence  suspension  theory  (O'Brien  1933,  Christiansen  1935,  Richardson 
1937,  and  Rouse  1937)  and  the  Prandtl-von  Karman  velocity  distribution 
equation  (Equation  2.2),  Lane  and  Kalinske  (1941)  applied  these  theories 
to  determine  suspended-sediment  load  using  a  point  suspended-sediment 
sample.  Lane  and  Kalinske  proposed  that  the  point  sample  could  be 
obtained  from  any  level  in  a  vertical  section;  however,  they  suggested 
that  the  sample  should  be  collected  from  the  lower  half-depth. 

Lane  and  Kalinske  applied  the  Prandtl-von  Karman  logarithmic  func¬ 
tion  for  the  vertical  velocity  distribution,  Equation  2.7.  Lane  and 
Kalinske  modified  Equation  2.7  using  von  Karman 's  constant  k  equal  to 
0.4  and  applying  Manning's  formula  for  hydraulic  roughness 


v  .  LdM  k2/3s1/2 

n 


2.39 


where  n  is  the  Manning's  roughness  coefficient.  The  resulting  form  of 
Equation  2.7  for  a  unit  width  of  a  wide  cross  section  (dw  R)  is 


Lane  and  Kallnske  substituted  Equations  2.40  and  2.24  (a  sediment 


distribution  equation)  Into  the  Integral 

d 

q  -  [  C  V  dy 

8  J  y  y 

0 

to  obtain 

r v 

q.  “  VCd  J  ^  exp 
8  7  0  V 

y 

where  q  Is  the  total  suspended-sediment  load  per  unit  width  and  C 
8  8 

Is  the  sampled  sediment  concentration  from  any  point  In  the  vertical. 
By  taking  exp  [lS(wa) /(UAd) ]  out  of  the  Integral,  Lane  and  Kallnske 
showed 


-  t5fe) 


(y  -  a) 


dy 


2.42 


q 


s 


qC^j  exp 


2.43 


where  q  Is  the  streamflow  In  cubic  feet  per  second  per  unit  width  and 
P^  is  a  function  of  w/U^  and  the  relative  roughness  n/d*^  .  Lane 
and  Kallnske  provided  a  graphical  solution  for  P^  . 

2. 2. 1.2  Brooks.  Since  the  research  of  Lane  and  Kallnske,  the 
advancement  of  determining  sediment  load  from  point  samples  of  suspended 
sediment  using  theoretically  based  concepts  received  little  attention 
until  the  research  efforts  of  Brooks  (1963).  Assuming  the  velocity 


distribution  in  Equation  2.7  and  the  concentration  distribution 
portrayed  in  Equations  2.26  and  2.27  (with  6*1),  Brooks  obtained 


where  is  the  aiddepth  sediment  concentration  at  y  *  d/2  and 

is  the  lower  limit  of  integration. 

Brooks  proposed  three  possible  ways  to  select  the  lower  limit  of 
integration,  y^  .  The  first  limit  is  Einstein's  (1950)  suggestion  that 
sediment  is  not  suspended  closer  to  the  bed  than  two  particle  diameters 
(2D).  Another  possible  limit  is  the  height  above  the  bed  where 
Vv  *  0  ,  according  to  the  logarithmic  velocity  profile.  This  second 
limit  prevents  unintentional  negative  contributions  to  the  suspended- 

sediment  load  Integral,  q  -  r  V  C  dy  ,  if  y_  is  greater  than  2D. 
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A  third  possible  limit  is  when  the  concentration,  determined  by  the 
extrapolation  of  the  suspended  load  equation,  is  equal  to  the  concentra¬ 
tion  of  the  size  fraction  in  the  bed.  Brooks  recommended  the  largest 
height  of  these  three  limits  as  the  most  reasonable  choice  for  the  limit 
to  integration. 

2.2.2  Conventional  Suspended- 
Sedlment  Measurement 

Following  the  development  and  improvement  of  depth-integrated 
suspended-sediment  samplers  in  the  1940' s,  the  convention  for  obtaining 
suspended-sediment  data  became  depth- Integrated  sampling.  Depth- 
integrated  sampling  is  limited  to  sampling  the  stream  section  referred 
to  as  the  measured  zone.  The  measured  zone,  under  most  conditions,  is 


within  the  region  of  suspended-sediment  transport.  The  unmeasured  zone 
contains  both  suspended  load  and  bed  load.  To  Include  the  unmeasured 
sediment  load  in  total  sediment  load  computations,  various  methods  have 
been  developed.  The  most  common  of  these  procedures  are  the  modified 
Einstein  procedure  (Colby  and  Hembree  1955) ,  the  Colby  method  (Colby 
1957),  and  the  Toffaletl  method  (Toffaleti  1969).  Of  these,  the  modi¬ 
fied  Einstein  procedure  has  received  the  most  attention  for  sediment 
load  measurements  (Schroeder  and  Hembree  1956,  Colby  and  Hubbell  1961, 
Lara  1966,  Bowie,  Bolton,  and  Murphree  1972,  Burkham  and  Dawdy  1980,  and 
Shen  and  Hung  1983) . 

2.2.2. 1  Modified  Einstein  Procedure.  The  modified  Einstein  pro¬ 
cedure  was  developed  by  Colby  and  Hembree  (1955)  to  determine  a  chan¬ 
nel's  total  sediment  load  through  the  combination  of  conventional 
suspended- sediment  discharge  measurements  and  a  bed  material  transport 
formula.  The  modified  Einstein  procedure  enables  the  computation  of  a 
channel's  unmeasured  load  which  can  be  added  to  a  measured  load  for  the 
estimation  of  total  load.  The  procedure  accomplishes  this  objective 
using  depth-integrated  suspended-sediment  samples  and  Einstein's  bed¬ 
load  function  (1950).  Data  required  to  use  the  procedure  include 
depth-integrated  suspended-sediment  samples,  streamflow  measurements, 
bed  material  gradations,  and  water  temperature.  Computational  compon¬ 
ents  of  the  modified  Einstein  procedure  include  Einstein's  modifications 
to  Keulegan's  open  channel  velocity  relationships  (Equations  2.18  and 
2.19  with  k  •  D._  and  flow  depth  replacing  hydraulic  radius).  House's 
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sediment  concentration  profile  relation  (Equation  2.26),  and  a  modified 
Einstein's  bed  load  function 
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,  **  .  3/2-3/2 f  ,,1/2 

"  r  Vsg  D  (ps-x) 


2.45 


Colby  and  Hembree  arbitrarily  divided  of  Equation  2.45  by  two  In 

order  for  their  sediment  load  calculations  to  agree  with  their  total 
sediment  load  measurements.  Substituting  Equations  2.18,  2.26,  and  2.45 
into  Equation  2.1  produces  the  sediment  load  equation  for  size  intervals 


Mt 


5.75U.  log 


M 


d?  +  Vb 


2.46 


where  1^  is  the  fraction  of  total  load  in  a  given  size  range  and  lfi 
is  the  fraction  of  bed  load  in  a  given  size  range.  Einstein  (1950) 
transformed  Equation  2.46  to 


*TqT 


11.6U.C  a 

*  a 


2.47 


in  which 


2.48 


2.49 


and  A  is  a  dimensionless  expression  equal  to  a/d  . 
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According  to  Einstein  (1950),  the  product  11. 6U*  is  the  effective 
flow  velocity  for  transporting  the  bed  material  in  the  bed  layer. 
Therefore,  in  addition  to  Equation  2.45,  bed  material  transport  in  the 
bed  layer  may  be  expressed  by 


Vb  -  ll*6U*Caa 


2.50 


Substitution  of  Equation  2.50  into  Equation  2.47  leads  to 


V>T  ‘  V»  (PIi  *h+l) 


2.51 


in  which 


P  -  2.303 


log 
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Before  Equation  2.51  can  be  applied  to  determine  the  total  sediment 
load,  the  value  for  z  must  be  determined.  The  value  for  z  is 
obtained  by  trial  and  use  of  the  measured  suspended-sediment  discharge 
for  the  dominant  size  interval  of  the  measured  zone. 

Colby  and  Hembree  reasoned  that  the  ratio  of  the  suspended-sediment 
load  of  the  measured  zone  to  the  suspended-sediment  load  for  the 

entire  suspended-sediment  load  zone  Qg  could  be  represented  by 


q;  pj{  + 

Qs  “  PJj  ♦  J2 


2.53 


where 


J 


2.54 


2.55 


and  a  single  prime  mark  designates  association  with  the  measured  zone. 
Equation  2.53  is  rewritten  with  the  substitution  of  the  relation  for  the 
total  suspended-sediment  load, 

Qs  -  WPI1  +  V  2-56 

where  i^Qg  is  the  sediment  discharge  through  the  bed  layer  of 
particles  of  a  given  size  range.  Substitution  of  Equation  2.56  into 
Equation  2.53  leads  to 


Qs  ll 

rr "  r(PJl +  JP 

B^B  1 


2.57 


where  Q'  is  determined  from  measurements  and  i_Q_  is  computed  using 
the  product  of  channel  width  and  Equation  2.45.  Equation  2.57  is  used 
to  compute  z  by  trial  for  the  dominant  size  interval  of  the  bed  load 
and  suspended  load.  The  values  of  z  for  the  other  size  ranges  are 
computed  using  the  following  relation  determined  by  Colby  and  Hembree: 
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z 


2.58 


where 

zy  *  z  for  size  ranges  other  than  the  reference  size 

z  *  z  for  the  reference  size 
r 

w^  *  w  for  size  ranges  other  than  the  reference  size 

w  ■  w  for  the  reference  size 
r 

2. 2. 2. 2  Colby's  Method.  Colby  (1957)  established  several  empiri¬ 
cal  relationships  for  estimating  total  sediment  discharge  using  the  same 
type  of  stream  measurements  used  in  developing  the  modified  Einstein 
procedure.  Colby's  procedure  is  simpler  to  apply  than  the  modified 
Einstein  procedure,  although  it  does  not  compute  the  sediment  load  for 
individual  size  fractions  and  it  is  sensitive  to  the  mean  velocity. 

Data  required  for  this  empirical  procedure  which  estimates  unmeasured 
sediment  load  are  the  mean  flow  velocity,  the  stream  width,  the  mean 
flow  depth,  and  the  measured  mean  suspended-sediment  discharge  concen¬ 
tration.  The  total  sediment  discharge  is  the  sum  of  the  measured  and 
unmeasured  sediment  loads. 

2. 2. 2. 3  Toffaleti's  Method.  Toffaleti  (1969)  suggested  that  his 
procedure  "offers  a  means  of  reducing  the  cost,  time,  and  manpower 
requirements. . .in  connection  with  continuing  sediment  observation  sta¬ 
tions."  He  further  suggested  that  his  procedure  could  be  used  to  com¬ 
pute  total  sand  load  which  can  be  added  to  wash  load  measured  from  three 
combined  samples  to  determine  total  sediment  load.  Wash  load  is  that 
part  of  the  total  sediment  discharge  that  is  composed  of  particle  sizes 
finer  than  those  represented  in  the  channel  bed. 
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Toffaleti  deviated  from  other  total  sand  load  computation  proce¬ 
dures  with  his  division  of  the  flow  depth  into  four  zones  of  sediment 
transport:  (1)  the  bed  zone  of  relative  thickness  2D/d  ;  (2)  the  lower 

suspended  zone  from  2D/d  to  y/d  ■  1/11 . 2 A  ;  (3)  the  middle  suspended 
zone  from  y/d  ■  1/11.24  to  y/d  ■  1/2.5  ;  and  (4)  the  upper  suspended 
zone  from  y/d  *  1/2.5  to  the  flow  surface.  For  each  of  the  three 
suspended  load  zones,  the  sediment  concentration  distribution  of  each 
sediment  size  fraction  is  expressed  as  power  relations.  Toffaleti' s 
relations  for  the  upper,  middle,  and  lower  zones,  respectively,  are 
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where 
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C  -  260.67  -  0.667T  2.63 

z 


25 


where 

Cu,Cm,Ct  ”  sediment  concentration  at  the  lower  boundaries  of  the 
upper,  middle,  and  lower  suspended  load  zones, 
respectively 

Sw  ■  water-surface  slope 
T  -  water  temperature 

also  equates  to  the  bed  layer  concentration. 

The  velocity  profile  for  the  flow  depth  is  represented  by  the  power 
relation 


+  s>*Ji) v 


where  the  exponent  z^  is  the  empirical  relation 


z  -  0.1198  +  0.00048T 
v 


Toffaleti  developed  an  empirical  relationship  for  the  sand  load  in 
the  lower  suspended  zone  to  make  his  sediment  load  computation  procedure 
complete.  This  empirical  relationship  provides  a  necessary  boundary 
condition  to  compute  the  sand  load  for  each  of  the  other  transport 
zones. 

The  following  input  is  required  for  the  Toffaleti  model:  the  aver¬ 
age  velocity  in  the  vertical,  the  hydraulic  radius,  water  temperature, 
stream  width,  bed  material  size  fraction,  Dgj  of  the  bed  material, 
water-surface  slope,  and  the  kinematic  viscosity  of  the  water-sediment 
mixture. 


2. 2. 2. 4  Burlchaa  and  Bawdy* s  Method.  Burkham  and  Dawdy  (1980) 

revised  Che  modified  Einstein  procedure  (Colby  and  Hembree  1955)  by 
eliminating  some  of  the  empirical  adjustments.  Burkham  and  Dawdy* s 
revision  eliminates  the  need  to  arbitrarily  divide  Einstein's  bed  load 
transport  intensity  factor,  of  Equation  2.28,  by  2  (Equation  2.45). 

Values  for  are  taken  directly  from  Einstein's  curve  relating 

to  (Einstein  1950,  Figure  10),  and  ^  is  determined  by  Equa¬ 

tion  2.34.  The  equivalent  sand  roughness  k  of  Equations  2.18,  2.19, 

s 

and  2.52  was  determined  to  be  5.5D^  rather  than  (Burkham  and 

Dawdy  1976).  Vlth  these  variations,  the  computation  steps  parallel 
those  of  Colby  and  Hembree  (1955). 

Burkham  and  Dawdy  compared  their  revision  of  the  modified  Einstein 
procedure  to  the  original  modified  Einstein  procedure  using  data  sets 
from  the  following  sites:  Niobrara  River  near  Cody,  Nebraska;  Middle 
Loup  River  at  Dunning,  Nebraska;  Fivemile  Creek  near  Riverton,  Wyoming; 
Fivemile  Creek  near  Shoshonl,  Wyoming;  and  Rio  Grande  conveyance  channel 
near  Bernardo,  New  Mexico.  For  these  data  sets,  Burkham  and  Dawdy  found 
no  significant  difference  in  the  sediment  discharges  computed  by  these 
two  procedures.  However,  both  procedures  did  significantly  under¬ 
estimate  sediment  loads  on  the  Rio  Grande  conveyance  channel  for  periods 
when  the  bed  form  was  flat  or  in  transition.  Also,  these  two  procedures 
significantly  overestimated  sediment  loads  for  two  sites  on  Fivemile 
Creek  in  Wyoming  where  Burkham  and  Dawdy  believe  the  bed  forms  were 
probably  dunes,  flat,  or  in  transition  at  different  times. 

2. 2. 2. 5  Shen  and  Hung's  Method.  Shen  and  Hung  (1983)  revised  the 
modified  Einstein  procedure  (Colby  and  Hembree  1955)  with  two  major 
modifications.  The  first  modification  applies  the  research  of  Lara 
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(1966)  which  proposed  a  nonlinear  relationship  between  z  and  w  . 

Shen  and  Hung  proposed  that  tht  relationship  between  z  and  w  can  be 
described  as 


a 

z  ■  djW 


2 


2.66 


where  and  a 2  are  regression  coefficients  determined  by  a  least 
squares  fit  using  field  data.  Equation  2.66  permits  a  variation  from 
Colby  and  Hembree's  proposed  relation  (Equation  2.58)  which  suggests 


and 


a2  -  0.7  2.68 

and  requires  the  selection  of  a  single  sediment  reference  size.  Shen 
and  Hung's  least  squares  fit,  Equation  2.66,  removes  the  subjectivity  in 
selecting  a  single  reference  size  and  permits  the  determination  of  a 
more  reliable  relationship.  Shen  and  Hung  show  that  z  does  not  vary 
with  w  raised  to  the  0.7  power.  The  other  modification  is  an  optimi¬ 
zation  procedure  to  minimize  the  difference  between  the  measured  and 
calculated  suspended-sediment  loads  in  the  sampled  zone. 

2. 2. 2. 6  Extrapolated  Data  Procedure.  Allen  and  Barnes  (1975) 
proposed  a  modification  to  conventional  suspended-sediment  load 


measurement.  Their  procedure  permits  the  direct  measurement  of  the 
zones  previously  classified  as  measured  and  unmeasured.  The  direct 
measurements  are  made  using  a  combination  of  depth-integrated  and  point 
sediment  samples  in  vertical  sections.  The  depth- Integrated  samples  are 
used  to  determine  the  sediment  load  In  the  measured  zone,  and  the  point 
samples  are  used  to  determine  the  suspended-sediment  load  in  the  zone 
not  measured  by  the  depth-integrated  sampler.  The  addition  of  these  two 
loads  equate  to  the  total  suspended-sediment  load.  Data  required  to 
determine  the  unmeasured  load  are  point  concentration  and  velocity  data 
from  two  or  more  depths  in  each  vertical  section.  These  point  concen¬ 
tration  and  velocity  data  are  plotted  on  log-log  plots  to  determine 
coefficients  for  the  concentration  and  velocity  relationships 


2.69 


and 


2.70 


respectively,  where 

Cq  -  concentration  at  (d  -  y)/y  ■  1.0 

z  •  log  Cj/Cq  where  Cj^  Is  the  concentration  at  (d  -  y)/y  ■  10.0 

Vq  -  velocity  at  y  *  1.0  on  the  distance  scale 

m  •  log  V,/V-  with  V.  being  the  velocity  at  y  -  10.0  on  the 

distance  scale 

The  product  of  Equations  2.69  and  2.70  Is  integrated  over  the  unmeasured 


zone  to  determine  the  unmeasured  load 
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The  lover  limit  of  integration  was  chosen  to  follow  the  suggestion  of 
Einstein  (1950)  and  the  upper  limit  of  integration  was  chosen  as  0.3  ft 
because  nearly  all  of  Allen  and  Barnes'  depth- integrated  samples  were 
taken  to  that  elevation. 

2.3  CHAPTER  SUMMARY 

The  review  of  the  literature  shows  that  the  vertical  velocity  dis¬ 
tribution,  the  sediment  concentration  distribution,  and  the  bed  load 
discharge  are  necessary  computational  components  for  determining  the 
total  sediment  load.  The  basis  for  determining  the  vertical  velocity 
distribution  by  most  of  Che  known  sediment  load  measurement  procedures 
is  the  Prandtl-von  Kerman  velocity  distribution  equation  (Equation  2.2). 
Most  of  the  reviewed  conventional  measurement  procedures  apply  the 
Keulegan  (1938)  form  of  the  vertical  velocity  distribution  equation  as 
it  was  modified  by  Einstein  (1950)  in  Equations  2.18  and  2.19.  The  con¬ 
centration  profile  relationship  applied  in  many  of  these  same  procedures 
is  Equation  2.26,  which  was  developed  by  Rouse  (1937).  The  bed  load 
discharge  Is  commonly  computed  with  a  form  of  the  Einstein  bed  load 
function  (Einstein  1950) ,  which  permits  the  determination  of  the  con¬ 
centration  at  the  lower  limit  of  the  suspended-sediment  zone. 


CHAPTER  3 


MODEL  FORMULATION 

3 . 1  NUMERICAL  COMPONENTS 

Conventional  total  sediment  load  measurement  procedures  use  depth- 
integrated  suspended-sediment  samples.  The  total  sediment  load  measure¬ 
ment  procedure  formulated  In  this  chapter  makes  use  of  point-source 
suspended-sediment  samples.  The  similarity  between  these  procedures  is 
the  use  of  relations  for  a  vertical  velocity  distribution,  a  suspended- 
sediment  concentration  profile,  and  bed  load  transport.  For  these  com¬ 
ponents,  the  proposed  procedure  applies  Einstein's  modifications  to 
Keulegan's  mean  and  point  velocity  relations  (Equations  2.18  and  2.19), 
Rouse's  sediment  concentration  distribution  relation  (Equation  2.26), 
and  Burkham  and  Dawdy's  modified  form  of  Einstein's  bed  load  function 
(Equations  2.28  and  2.34).  The  primary  feature  that  sets  the  proposed 
procedure  apart  from  other  methods  is  the  direct  computation  of  z  for 
Equation  2.26  using  a  point  sample  of  suspended  sediment.  When  Equa¬ 
tion  2.26  is  applied,  z  is  the  only  unknown,  provided  the  suspended- 
sediment  concentration  is  known  at  a  known  depth  and  the  sediment  con¬ 
centration  at  the  lower  suspended-sediment  boundary  is  calculated. 
Section  3.1.2  and  Appendix  A  provide  greater  detail  on  the  calculation 
of  z  . 

3.1.1  Velocity  Distribution 

Function _ 

The  velocity  distribution  and  mean  velocity  equations  used  by  the 
reviewed  conventional  sediment  load  measurement  procedures  are 
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Equations  2.18  and  2.19,  respectively.  Variations  in  the  application  of 

these  two  equations  come  from  the  use  of  different  relations  for  k  . 

8 

Einstein  (1950)  suggested  k  -  D, c  as  a  good  rule  of  thumb  for 

s  o  J 

application  in  Equations  2.18  and  2.19.  Einstein's  rule  has  been 
applied  by  each  of  the  conventional  sediment  measurement  procedures 
using  these  equations  except  for  the  modification  proposed  by  Burkham 
and  Dawdy  (1980).  Burkham  and  Dawdy  (1976)  analyzed  the  data  of  Barnes 
(1967)  and  Limerlnos  (1970)  and  determined  that  k  -  5.5D,C  .  Burkham 

8  03 

and  Dawdy  (1976)  also  indicated  that  the  use  of  a  larger  referenced 

particle  size  improves  the  prediction  of  the  mean  velocity  and  velocity 

distribution.  If  a  larger  particle  size  is  to  be  used,  Burkham  and 

Dawdy  agree  with  the  use  of  kg  -  3.3Dg^  as  aPPHed  by  Limerlnos.  The 

Limerlnos  relationship  for  k  is  compatible  with  van  Ri_1n's  (1982) 

8 

review  of  published  k  relationships  and  recommendation  that 

8 

kg  -  3°9o  for  grain  roughness. 

Literature  shows  for  movable-bed  streams  that  k  is  variable  and 

s 

dependent  upon  flow  regime  (Nordin  and  Dempster  1963,  van  Rijn  1982). 

The  lower  limit  of  kg  results  from  grain  roughness  during  upper  regime 
flow  (i.e.,  plane  bed  and  antidunes),  and  the  upper  limit  of  kg  re¬ 
sults  from  both  grain  and  form  roughness  during  lower  regime  flow  (i.e., 
ripples  and  dunes) .  The  Importance  of  form  roughness  was  demonstrated 
by  Nordin  and  Dempster  (1963)  with  their  study  of  the  Middle  Rio  Grande, 
New  Mexico.  Nordin  and  Dempster  calculated  for  the  Middle  Rio  Grande  an 

average  maximum  k  of  400D..  . 

S  03 

The  velocity  distribution  of  the  proposed  procedure  uses  Equa¬ 
tions  2,18  and  2.19  as  in  most  of  the  conventional  sediment  measurement 
procedures  (Colby  and  Hembree  1955,  Burkham  and  Dawdy  1980,  and  Shen  and 
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Hung  1983) .  A  variation  in  the  proposed  procedure  from  the  conventional 
procedures  Is  the  provision  for  a  variable  kg  and  the  use  of  Dg^  of 
the  bed  material  as  a  reference  size.  A  variable  kg  makes  it  possible 
to  account  for  both  grain  and  form  roughness.  The  larger  referenced  bed 
material  size  Improves  the  prediction  of  the  mean  velocity  and  velocity 
distribution  as  indicated  by  Burkham  and  Dawdy  (1976).  The  reference 
size  multiplier  may  be  selected  based  on  the  judgment  of  the  investi¬ 
gator  or  determined  using  the  measured  stream  hydraulics  and  computa¬ 
tional  procedures  such  as  Einstein  (1950)  and  Brownlie  (1981,  1983). 

The  Einstein  and  Brownlie  procedures  predict  the  resistance  to  flow  over 
movable  boundaries  with  bed  forms.  The  use  of  Einstein's  and  Brownlie's 

procedures  to  determine  k  is  presented  in  Appendix  A. 

8 

3.1.2  Concentration  Profile 

The  suspended-sediment  concentration  profile  in  the  proposed  pro¬ 
cedure  is  determined  according  to  Equation  2.26.  The  use  of  Equa¬ 
tion  2.26  for  determining  the  concentration  profile  is  similar  to  most 
conventional  sediment  measurement  procedures  (Colby  and  Hembree  1955, 
Burkham  and  Dawdy  1980,  and  Shen  and  Hung  1983)  and  varies  only  in  the 
determination  of  z  .  The  calculation  for  z  is  based  upon  the  sampled 
suspended-sediment  concentration  and  the  calculation  of  the  bed  load 
concentration  using  Burkham  and  Dawdy' s  (1980)  modification  of 
Einstein's  bed  load  function.  Values  for  z  are  directly  computed 
using  the  modified  form  of  Equation  2.26 

/  \ 


where 


Cb  •  the  bed  load  sediment  concentration, 

C  •  the  concentration  of  the  suspended  sediment  at  distance  y 
"  above  the  channel  bed, 

yb  ■  the  height  above  the  channel  bed  for  the  lower  boundary  of 
suspended  sediment 

y  »  the  height  of  suspended-sediment  point  sample  above  the 
^  channel  bed 


The  lower  boundary  for  suspended  sediment  is  set  at  the  largest  of 
two  possible  limits.  The  first  limit  is  Einstein's  (1950)  suggestion 
that  sediment  is  not  suspended  closer  to  the  bed  than  two  particle  diam¬ 
eters  (2D).  The  ocher  limit  is  the  height  above  the  bed  where  V  »  0  , 
according  to  Equation  2.18.  This  latter  limit  prevents  unintentional 
negative  contributions  to  the  suspended-sediment  load  Integral, 

q  »  |  VC  dy  ,  if  y,  is  greater  than  2D. 

3  Jv  y  y  b 

yb 


3.1.3  Bed  Load  Function 

Burkham  and  Dawdy's  (1980)  modified  form  of  Einstein's  bed  load 
function  (Equations  2.28  and  2.34)  is  applied  to  compute  the  bed  load. 
Burkham  and  Dawdy's  modification  is  chosen  because  it  permits  the  com¬ 
putation  of  bed  load  for  individual  sediment  size  fractions  and  assists 
in  the  elimination  of  arbitrarily  dividing  by  2  as  in  Equation  2.45 

(compare  to  Einstein's  equation.  Equation  2.28). 

To  apply  Equation  2.34,  the  product  of  hydraulic  radius  and  energy 
slope  RS  must  be  determined.  Einstein  (1950)  recommended  that  the 
product  RS  should  be  separated  into  R'S  and  R"S,  where  R  ■  R'  +  R"  , 
and  R"  is  the  hydraulic  radius  due  to  form  roughness.  Einstein  based 
his  recommendation  for  the  separation  of  R  into  two  parts  on  reasoning 
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that  flow  energy  is  transformed  Into  energy  of  turbulence  which  may  be 
separated  into  two  forms.  The  first  form  of  turbulence  is  caused  by  the 
rough  wall  and  stays  in  the  vicinity  of  the  boundary  material.  The 
second  form  of  turbulence  corresponds  to  shape  resistance  and  occurs  at 
a  considerable  distance  away  from  the  grains.  Einstein  further  reasoned 
that  the  energy  due  to  the  first  form  affects  the  bed  load  motion  and 
the  energy  due  to  the  second  form  does  not  contribute  to  the  bed  load 
motion. 

The  division  of  R  into  R'  and  R"  leads  to  added  numerical 
calculations  and  was  not  applied  by  Colby  and  Hembree  (1955)  nor  any  of 
the  subsequent  sediment  measurement  methods.  The  proposed  procedure 
applies  Einstein's  reasoning  to  separate  the  flow  energy;  however,  the 
product  RS  ,  is  separated  into  RS'  and  RS”  (as  suggested  by  Meyer- 
Peter  and  Muller  1948  and  Vanonl  and  Brooks  1957) ,  where 
S  -  s'  +  S" 

S'  •  the  energy  slope  due  to  grain  roughness 

S"  *  the  energy  slope  due  to  form  roughness 
Equation  2.34  may  now  be  written  as 


K  -  1.65(51)*,  \T^J  r§i 


and  the  product  RS'  is  computed  using  the  modified  form  of 


Equation  2.19, 


=  -  5.75  log  ^12.27 


where  k'  is  the  effective  roughness  element  height  due  to  grain  rough- 
s 

ness  and  x*  is  a  function  of  k'/6'  which  is  presented  in  Figure  A  of 

s 

Einstein  (1950),  and  5'  ■  (11.6v)/U^  .  Note  that  Equation  3.3  provides 
the  same  computational  results  as  the  conventional  sediment  measurement 
procedures  when  bed  forms  are  insignificant. 

3 . 2  MODEL  DESCRIPTION 

The  computational  relations  in  the  proposed  procedure  applv 
Einstein's  modifications  to  Keulegan's  (1938)  mean  and  point  velocity 
relations  (Equations  2.18,  2.19,  and  3.3),  the  suspended-sediment  con¬ 
centration  distribution  Equation  2.26,  and  Burkham  and  Dawdy's  (1980) 
modified  form  of  Einstein's  bed  load  function  (Equations  2.28  and  3.2). 
The  basic  data  requirements,  computation  steps,  and  computation  results 
are  similar  to  the  modified  Einstein  procedure  and  its  revisions.  The 
major  difference  in  the  proposed  procedure  is  the  collection  of 
suspended-sediment  data  from  a  point  sample  rather  than  from  a  depth- 
integrated  sample.  The  data  requirements  and  computation  steps  for  the 
proposed  procedure  are  summarized  in  the  following  sections  of  this 
chapter. 

3.2.1  Basic  Data 

To  effectively  apply  the  proposed  procedure,  the  basic  data  for 
each  total  sediment  load  measurement  are  to  represent  the  field 
conditions  at  the  time  of  suspended-sediment  sampling.  The  basic  data 
for  the  proposed  procedure  are  the  streamflow  for  the  sampled  cross 
section;  the  sampled  cross  section  width,  cross-sectional  area,  and 
hydraulic  radius;  the  flow  depth  of  the  sampled  vertical;  the  height 
above  the  channel  bed  where  the  point  sample  was  obtained;  the  con¬ 
centration  of  the  suspended-sediment  sample;  the  fraction  of  bed 


material  in  individual  size  intervals;  the  fraction  of  suspended  mate¬ 
rial  in  individual  size  Intervals;  and  the  water  temperature. 

3.2.2  Computation  Steps 

The  general  application  of  the  computation  components  and  field 
data  for  determining  total  sediment  load  is  presented  in  detail  in 
Appendix  A.  The  Fortram  77  computer  code  for  these  computations  is 
provided  In  Appendix  B.  A  summary  of  the  computations  used  with  the 
basic  data  for  the  proposed  procedure  is  presented  in  step  form  as 
follows: 

Step  1.  Compute  the  bed  load  for  a  unit  width  of  flow  for  each  bed 
material  size  interval  (Equations  2.28  and  2.34-2.38). 

Step  2.  Using  the  results  of  step  1,  compute  the  concentration  in  the 
bed  load  zone  for  each  bed  material  size  interval 
(Equation  A. 9). 

Step  3.  Determine  the  exponent  z  of  Equation  2.26  for  each  sediment 

size  interval  using  the  concentration  of  the  suspended-sediment 
sample  and  the  concentration  in  the  bed  load  zone  determined 
in  step  2. 

Step  4.  Using  the  results  of  step  3,  develop  a  linear  regressed  rela¬ 
tionship  for  z  as  a  function  of  particle  fall  velocity  (Equa¬ 
tion  A. 12) . 

Step  5.  Recompute  the  bed  load  zone  concentration  for  each  sediment 

size  interval  using  the  sampled  suspended-sediment  concentra¬ 
tions,  the  relationship  for  z  developed  in  step  4,  and  Equa- 


Step  6.  Recompute  the  bed  load  for  a  unit  width  of  flow  for  each 

sampled  suspended-sediment  size  interval  using  the  results  of 
step  5  (Equation  A. 15). 

Step  7.  Determine  the  total  sediment  load  for  a  unit  width  of  flow  for 
each  sediment  size  interval  using  Equation  2.51,  the  computed 
bed  load,  and  the  relationship  for  z  developed  in  step  4. 

Step  8.  Compute  the  total  sediment  load  for  the  channel  by  multiplying 
the  sum  of  the  results  from  step  7  by  the  width  of  the  sampled 
cross  section  (Equation  A. 18). 

3.3  CHAPTER  SUMMARY 

A  procedure  for  measuring  total  sediment  load  using  suspended- 
sediment  data  collected  at  one  point  has  been  developed.  The  procedure 
advances  the  research  efforts  of  Lane  and  Kallnske  (1941)  and  Brooks 
(1963),  eliminates  the  need  for  a  correlation  between  point-source  and 
depth- integrated  suspended-sediment  samples,  and  permits  the  direct  com¬ 
putation  of  total  sediment  load  using  point-source  suspended-sediment 
data.  The  proposed  procedure  applies  Einstein's  (1950)  modification  to 
Keulegan's  (1938)  mean  and  point  velocity  relations,  the  suspended- 
sediment  concentration  distribution  equation  developed  by  Rouse  (1937), 
and  Burkham  and  Dawdy's  (1980)  modified  form  of  Einstein's  bed  load 
function.  The  data  required  for  application  of  the  proposed  procedure 
include  descriptors  of  the  streamflow,  the  channel  geometry,  the  bed 
material,  and  the  point-source  suspended-sediment  sample. 


CHAPTER  4 


MODEL  EVALUATION 

4 . 1  GENERAL 

This  chapter  compares  the  proposed  sediment  load  measurement  proce¬ 
dure  to  the  total  sediment  load  field  data  and  the  modified  Einstein 
procedure.  The  computational  step  solution  and  computer  code  listing 
for  the  proposed  procedure  are  provided  in  Appendices  A  and  B,  respec¬ 
tively.  The  proposed  procedure  and  computer  code  are  identified  as  TSL 
(Total  Sediment  Load) . 

4.2  BASIC  DATA 

The  TSL  procedure  is  compared  to  observed  total  sediment  load  data 
and  the  modified  Einstein  procedure  using  data  from  four  selected  river 
reaches:  the  Niobrara  River  at  Cody,  Nebraska  (Colby  and  Hembree  1955); 
Fivemile  Creek  at  Shoshoni,  Wyoming  (Colby,  Hembree,  and  Rainwater 
1956);  the  Middle  Loup  River  near  Dunning,  Nebraska  (Hubbell  and 
Matejka  1959);  and  the  Rio  Grande  conveyance  channel  near  Bernardo, 

New  Mexico  (Culbertson,  Scott,  and  Bennett  1972).  These  data  sets 
were  selected  bee  ruse  of  their  use  in  the  development  and  evaluation  of 
the  modified  Einstein  procedure.  Also,  the  selected  data  sets  provide 
the  suspended-sediment  point  samples  required  for  application  of  the  TSL 
procedure. 

4.3  COMPARISON  TO  TOTAL 
SEDIMENT  LOAD  MEASUREMENTS 

The  referenced  data  sets  provide  sediment  data  for  both  normal  and 
total  sediment  load  sampling  stations.  The  "normal"  sediment  sampling 
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stations  are  located  in  channel  reaches  with  a  uniform  cross  section  and 
uniform  flow.  The  total  sediment  load  sampling  stations  are  located 
downstream  of  the  normal  sampling  stations  at  natural  channel  constric¬ 
tions  or  man-made  turbulence  flumes  wherein  depth-integrated  suspended- 
sediment  samples  are  assumed  to  represent  all  of  the  sediment  in 
transport . 

4.3.1  Niobrara  River 

Colby  and  Hembree  (1955)  describe  the  study  reaches,  the  data  col¬ 
lection,  and  the  data  for  the  Niobrara  River  near  Cody,  Nebraska.  Colby 
and  Hembree  used  the  Niobrara  River  data  to  develop  the  modified 
Einstein  procedure.  The  data  Include  depth- Integrated  suspended- 
sediment  samples  which  were  concurrently  collected  at  the  normal  and 
contracted  sections  for  comparison  of  the  modified  Einstein  procedure  to 
the  total  sediment  load.  The  data  also  include  point-integrated 
suspended-sediment  samples  which  were  collected  to  assess  the  suspended- 
sediment  concentration  distribution.  The  point  samples,  however,  were 
not  collected  concurrently  with  the  total  sediment  load  measurements. 
Therefore,  for  this  study,  comparisons  of  the  TSL  procedure  to  the  total 
sediment  load  measurements  are  made  by  comparing  sediment  load  versus 
streamflow  relationships. 

Total  sediment  loads  observed  at  the  contracted  section  of  the 
Niobrara  River  are  plotted  against  streamflow  in  Figure  4.1.  The 
scatter  of  points  in  Figure  4.1  shows  a  trend  of  Increasing  sediment 
load  with  increasing  streamflow  and  can  be  regressed  to  the  function, 
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Figure  4.1.  Observed  total  sediment  load  versus  streamflov 
for  the  Niobrara  River.  The  curve  represents 
Equation  4.1. 


4.1 


Qj  -  0.003672Q2  +  6.639Q  -  1282.0 

where  Is  the  total  sediment  load  In  tons  per  day  and  Q  is  the 
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streamflov  in  cubic  feet  per  second.  The  coefficient  of  correlation  r 
for  Equation  4.1  is  0.89. 

Comparison  of  the  observed  total  sediment  loads  and  the  total  sedi¬ 
ment  loads  estimated  by  Equation  4.1  is  provided  in  Figure  4.2.  Fig¬ 
ure  4.2  shows  87  percent  of  the  observed  sediment  loads  to  be  within 
150  percent  of  the  sediment  load  estimated  by  Equation  4.1.  This  per¬ 
centage  difference  is  an  Indicator  of  the  accuracy  of  an  attempt  to 
determine  a  sediment  load  relationship  for  the  Niobrara  River. 

Total  sediment  loads  computed  by  the  TSL  procedure  using  point 
samples  at  the  normal  section  are  plotted  against  streamflow  in  Fig¬ 
ure  4.3.  The  scatter  of  points  in  Figure  4.3  shows  a  trend  of  increas¬ 
ing  sediment  load  with  Increasing  streamflow.  The  distribution  of  sedi¬ 
ment  load  magnitudes  for  a  given  range  of  streamflow  in  Figure  4.3  is 
greater  than  the  distribution  in  Figure  4.1.  This  increased  distribu¬ 
tion  may  result  from  more  samples  or  the  sensitivity  of  the  point 
sampler' 8  location.  The  point  samples  for  Figure  4.3  were  taken  from 
the  lower  14  to  90  percent  of  the  flow  depths  (0.5  to  3.1  ft  above  the 
channel  bed) . 

Lane  and  Kallnske  (1941)  recommended  that  suspended-sediment  point 
samples  should  be  collected  from  the  lower  half-depth.  Point  sample 
measurements  used  to  further  assess  the  TSL  procedure  were  from  the 
lower  third  flow  depth  to  be  well  within  Lane  and  Kalinske's  recommenda¬ 
tion.  For  example,  data  for  samples  collected  in  the  lower  third  of  the 
flow  depths  (0.5  ft  above  the  channel  bed)  are  plotted  in  Figure  4.4. 
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Figure  A. 2.  Observed  sediment  load  versus  estimated  sediment  . 
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Figure  4.3.  Total  sediment  load  determined  by  TSL  versus  stresmflow 
for  the  Niobrara  River  using  point  samples  from  the 
lower  14-90  percent  flow  depth  end  the  Llmerlnos 
procedure  to  determine  k  .  The  curve  represents 
Equation  4.1.  9 
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.4.  Total  sediment  load  determined  by  TSL  versus  streamflov 
for  the  Niobrara  River  using  point  samples  from  the 
lover  one-third  flow  depth  and  the  Llmerlnos  procedure 
to  determine  k  .  The  curve  represents  Equation  4.1. 


The  scatter  of  points  in  Figure  4.4  is  reduced  from  the  level  of  scatter 
in  Figure  4.3.  The  reduction  of  data  scatter  in  Figure  4.4  may  be  at¬ 
tributed  to  greater  reliability  in  sampling  sediment  concentrations  in 
the  lower  flow  depth  as  Lane  and  Kallnske  recommended.  At  locations 
further  from  the  channel  bed  the  sediment  concentration  and  flow 
velocity  may  deviate  from  the  theoretical  profiles  due  to  factors  such 
as  macroturbulence. 

The  data  of  Figure  4.4  are  replotted  in  Figure  4.5  to  assess  the 
relationship  between  the  TSL  procedure  and  the  estimated  total  sediment 
loads  as  represented  by  Equation  4.1.  Figure  4.5  shows  73  percent  of 
the  sediment  loads  determined  by  TSL  to  be  within  ±50  percent  of  the 
sediment  load  estimated  by  Equation  4.1.  These  values  are  a  measure  of 
the  accuracy  of  TSL  only  if  Equation  4.1  is  the  true  sediment  load  re¬ 
lationship  for  the  Niobrara  River  near  Cody,  Nebraska. 

Figures  4.4  and  4.5  show  that  9  percent  of  the  sediment  loads 
determined  by  the  TSL  procedure  are  at  least  100  percent  greater  than 
the  observed  sediment  loads  as  represented  by  Equation  4.1.  The  largest 
of  these  deviations  occur  at  the  larger  streamflows.  The  assessment  of 
these  deviations  is  limited  due  to  the  limited  observed  data  in  the 
larger  streamflow  range.  If  additional  observed  data  are  acquired  for 
the  larger  streamflows,  the  deviations  of  the  TSL  sediment  loads  from  a 
corrected  representative  sediment  load  relationship  may  not  be  as  sig¬ 
nificant  as  they  now  appear.  For  example,  these  deviations  would  be 
smaller  if  the  corrected  sediment  load  relationship  would  indicate 
larger  sediment  loads  for  the  larger  streamflows.  Also,  if  additional 
observed  data  would  display  greater  variability  about  a  regressed 
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Figure  4.5.  Sediment  load  determined  by  TSL  using  point  samples 

from  the  lower  one-third  flow  depth  and  the  Llmerlnos 
procedure  to  determine  k  versus  estimated  sediment 
load  for  the  Niobrara  Rlvlr. 


sediment  load  function,  the  significance  of  the  variability  In  the  TSL 


sediment  loads  would  be  lessened. 


Assuming  the  TSL  deviations  are  departures  from  the  true  sediment 


load  relationship  for  the  Niobrara  River,  the  deviations  may  be  attri¬ 


buted  to  measurement  error  or  spikes  In  the  instantaneous  point  concen¬ 


trations.  Measurement  errors  or  Instantaneous  spikes  would  be  removed 


when  TSL  measurements  are  collected  with  near-continuous  sampling.  Con¬ 


tinuous  point-source  suspended-sediment  sampling  and  the  TSL  procedure 


.  \ 


enable  tine  Integration  for  the  measurement  of  the  true  total  sediment 
load.  Time  Integration  results  in  an  average  of  consecutive  total  sedi¬ 
ment  load  measurements  about  a  near-constant  streamflov. 

Figure  4.6  provides  Insight  into  the  results  of  time  integration 
for  the  determination  of  a  river's  total  sediment  load.  The  sediment 
loads  of  Figure  4.6  are  average  TSL  sediment  loads  for  100-cfs- 
streamflow  class  intervals.  The  sediment  loads  shown  in  Figure  4.6  are 
within  the  ±50  percent  deviation  of  the  observed  sediment  loads  about 
Equation  4.1  as  shown  in  Figure  4.2.  Therefore ,  the  deviations  of 
Figures  4.4,  4.5,  and  4.6  are  acceptable. 

The  Niobrara  River  data  were  used  to  assess  the  Importance  of 

estimating  the  roughness  element  height  k  due  to  both  grain  and  form 
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roughness.  Section  3.1.1  discussed  the  Importance  of  estimating  the 
value  for  kg  to  be  used  in  the  velocity  distribution  function.  Equa¬ 
tion  2.18.  The  Llmerlnos  (1970),  the  Einstein  (1950),  and  the  Brownlie 
(1981,  1983)  roughness  prediction  procedures  may  be  used  to  determine 
kg  and  are  provided  as  options  in  the  computer  program,  TSL.  The  sedi¬ 
ment  loads  determined  by  TSL  with  these  procedures  are  presented  in  Fig¬ 
ures  4.4,  4.5,  4.7,  4.8,  4.9,  and  4.10.  These  figures  show  similar 
degrees  of  scatter  in  the  higher  discharge  range  about  the  observed 
sediment  loads  and  those  calculated  by  Equation  4.1.  The  discussion 
relevant  to  the  scatter  in  the  higher  discharge  range  remains  the  same 
as  the  previous  discussion  on  Figures  4.4  and  4.5.  The  sediment  loads 
in  the  lower  discharge  range  are  affected  more  by  the  differing  rough¬ 
ness  prediction  procedures  than  those  in  the  upper  discharge  range.  The 
different  sediment  load  estimates  in  the  lower  discharge  range  are 
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Figure  4.7. 


Total  sediment  load  determined  by  TSL  versus  streamflov 
for  the  Niobrara  River  using  point  samples  from  the 
lower  one-third  flow  depth  and  the  Brownlie  procedure 


to  determine 


The  curve  represents  Equation  4.1. 
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Sediment  load  determined  by  TSL  using  point  samples 
from  the  lower  one-third  flow  depth  and  the  Brownlie 
procedure  to  determine  k  versus  estimated  sediment 
load  for  the  Niobrara  Rlvfr. 
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Sediment  load  determined  by  TSL  using  point  samples 
from  the  lower  one-third  flow  depth  and  the  Einstein 
procedure  to  determine  k  versus  estimated  sediment 
load  for  the  Niobrara  Rivlr. 


attributed  to  the  Einstein  and  the  Brownlie  procedures  estimating  larger 
roughness  element  heights  than  the  Limerinos  procedure.  A  larger  rough¬ 
ness  element  height  prediction  leads  to  a  decreased  estimate  of  total 
sediment  load  for  the  same  point-source  suspended-sediment  sample.  In 
summary.  Figures  4.4,  4.5,  and  4.7-4.10  indicate  that  total  sediment 
load  measurement  by  TSL  is  sensitive  to  the  selection  of  the  roughness 
height  prediction  method.  The  analyses  of  these  figures  encourage  the 
selection  of  the  Limerinos  procedure. 

4.3.2  Fivemile  Creek. 

Colby,  Hembree,  and  Rainwater  (1956)  describe  the  study  reaches  and 
the  sediment  data  for  Fivemile  Creek  near  Shoshoni,  Wyoming.  Channel 
geometry  and  streamflow  rating  data  are  available  in  the  files  of  the 
US  Geological  Survey,  Cheyenne,  Wyoming.  Channel  geometry  and  stream- 
flow  rating  data  for  Fivemile  Creek  are  presented  in  Appendix  C  as  they 
were  used  for  this  analysis. 

Colby,  Hembree,  and  Rainwater's  report  presents  depth-  and  point- 
integrated  suspended-sediment  data  for  normal  and  contracted  sections. 
The  depth- integrated  suspended-sediment  data  at  the  contracted  section 
are  assumed  to  represent  all  of  the  sediment  in  transport.  The  point- 
integrated  suspended-sediment  data  were  not  collected  concurrently  with 
the  total  sediment  load  measurements.  Therefore,  comparison  of  the  TSL 
procedure  to  the  total  sediment  load  measurements  is  made  with  sediment 
load  versus  streamflow  relationships. 

Total  sediment  loads  observed  at  the  contracted  section  of  Fivemile 
Creek  are  plotted  against  streamflow  in  Figure  4.11.  The  scatter  of 
points  in  Figure  4.11  shows  a  trend  of  increasing  sediment  load  with 
increasing  streamflow.  A  single  function  for  the  plotted  data  is  not 


evident;  therefore,  results  of  the  TSL  procedure  are  compared  by  plot¬ 
ting  both  the  observed  and  TSL  data  sets  in  Figure  4.12.  The  suspended- 
sediment  point  samples  for  Figure  4.12  were  from  the  lower  10  to  33  per¬ 
cent  of  the  flow  depths  (0.2  to  0.4  ft  above  the  channel  bed).  Most  of 
the  sediment  load  data  in  Figure  4.12  indicate  good  agreement  between 
the  observed  and  TSL  data  sets.  Some  of  the  TSL  measurements  are  lower 
than  the  observed  sediment  loads  in  the  larger  streamflow  range.  These 
lower  sediment  load  measurements  may  represent  true  sediment  loads  at 
the  time  each  point-source  suspended-sediment  sample  was  obtained. 
Suspended-sediment  samples  for  the  observed  and  TSL  measurements  were 
not  concurrently  collected;  therefore,  this  latter  explanation  cannot  be 
further  assessed.  Other  possible  explanations  for  the  low  TSL  measure¬ 
ments  include  sampling  error  or  an  improper  calculation  by  TSL  for  bed 
load,  suspended-sediment  concentration  distribution,  or  velocity  pro¬ 
file.  The  possibility  of  improper  calculation  by  TSL  is  minimized  when 
consideration  is  given  to  the  other  TSL  data,  which  match  the  distribu¬ 
tion  of  the  observed  sediment  loads.  In  summary,  Figure  4.12  indicates 
good  agreement  between  the  observed  and  TSL  data  sets. 

4.3.3  Middle  Loup  River 

Hubbell  and  Matejka  (1959)  describe  the  study  reaches,  the  data 
collection,  and  the  data  for  the  Middle  Loup  River  at  Dunning,  Nebraska. 
Hubbell  and  Matejka's  report  provides  depth-  and  point-integrated 
suspended-sediment  data  for  normal  sections  and  a  turbulence  flume.  The 
depth- integrated  suspended-sediment  data  from  the  turbulence  flume  are 
assumed  to  represent  all  of  the  sediment  in  transport.  The  point- 
integrated  suspended-sediment  data  were  not  collected  concurrently  with 
the  total  sediment  load  measurements.  Therefore,  comparison  of  the  TSL 
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Figure  4.12.  Total  sediment  load  determined  by  TSL  versus  streamflov 
for  Flvemlle  Creek  using  point  samples  from  the  lover 
one-third  flow  depth  and  the  Limerlnos  procedure  to 
determine  k  . 


procedure  to  the  total  sediment  load  measurements  is  made  by  comparing 
sediment  load  versus  streamflov  relationships. 

Total  sediment  loads  observed  at  the  turbulence  flume  of  the  Middle 
Loup  River  are  plotted  against  streamflow  in  Figure  4.13.  Results  of 
the  TSL  procedure  for  the  reported  point-source  suspended-sediment  data 
are  assessed  by  plotting  the  observed  and  TSL  sediment  load  sets  in 
Figure  4.14.  The  suspended-sediment  point  samples  for  Figure  4.14  were 
from  the  lower  19  to  33  percent  of  the  flow  depths  (0.2  to  0.4  ft  above 
the  channel  bed).  Figure  4.14  shows  good  agreement  between  the  avail¬ 
able  observed  and  TSL  data  sets. 

4.3.4  Rio  Grande  Conveyance  Channel 

0 

Culbertson,  Scott,  and  Bennett  (1972)  describe  the  study  reaches, 
the  data  collection,  and  the  data  for  the  Rio  Grande  conveyance  channel 
near  Bernardo,  New  Mexico.  Culbertson,  Scott,  and  Bennett’s  report  pro¬ 
vides  depth-integrated  suspended-sediment  data  for  normal  sections  and  a 
weir  structure.  Point-integrated  suspended-sediment  data  for  a  narrow 
discharge  range  are  also  presented  for  the  sampled  sections  in  the 
report.  The  point  suspended-sediment  samples  were  not  collected  con¬ 
currently  with  the  total  sediment  load  measurements. 

Total  sediment  loads  observed  at  the  weir  structure  of  the  Rio 
Grande  conveyance  channel  are  plotted  against  streamflow  in  Figure  4.15. 
Results  of  the  TSL  procedure  for  suspended-sediment  point  samples  are 
compared  to  the  observed  data  in  Figure  4.16.  The  point  samples  for 
these  TSL  measurements  came  from  the  lower  6  to  32  percent  of  the  flow 
depths  (0.3  to  1.2  ft  above  the  channel  bed).  The  distribution  of  the 
TSL  sediment  loads  is  greater  in  Figure  4.16  than  the  observed  sediment 
loads.  The  greater  distribution  may  be  attributed  to  the  use  of 
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Figure  4.14.  Total  sediment  load  determined  by  TSL  versus  streamflow 
for  the  Middle  Loup  River  using  point  samples  from  the 
lower  one-third  flow  depth  at  section  C  and  the  Limerinos 
procedure  to  determine  k 
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point-source  suspended-sediment  samples  which  were  collected  too  close 
to  the  channel  bed.  Point  samples  collected  near  the  channel  bed  may 
Inadvertently  Include  bed  load  and  bed  material.  Also,  point  samples 
collected  close  to  the  channel  bed  may  exhibit  greater  variance  than 
other  point  samples  due  to  the  same  processes  that  influence  the 
stochastic  nature  of  bed  load  transport. 

The  use  of  point  samples  collected  near  the  channel  bed  was 
assessed  by  plotting  TSL  results  for  point  samples  taken  from  the  lower 
20  to  32  percent  of  the  flow  depths  (0.8  to  1.2  ft  above  the  channel 
bed) .  The  lower  limit  of  20  percent  was  selected  to  be  well  above  the 
previous  lower  limit  of  6  percent.  The  plot  of  TSL  measurements  using 
these  point  samples  is  presented  in  Figure  4.17.  Figure  4.17  shows 
better  agreement  between  the  observed  and  TSL  data  sets  than 
Figure  4.16. 

In  summary,  the  TSL  measurement  procedure  is  improved  when  the 
point  samples  are  obtained  in  the  lower  flow  depth  but  not  too  close  to 
the  channel  bed.  Point-source  suspended-sediment  samples  collected  in 
the  lower  20  to  32  percent  flow  depth  provide  reliable  data  for  the  TSL 
procedure. 

4.4  COMPARISON  TO  MODIFIED 

EINSTEIN  PROCEDURE  COMPUTATIONS 

One  of  the  objectives  of  this  research  was  to  compare  the  TSL  pro¬ 
cedure  to  the  modified  Einstein  procedure  (MEP) .  For  this  comparison, 
results  of  the  MEP  as  applied  to  the  published  data  are  plotted  with  the 
observed  total  sediment  load  measurements  in  Figures  4.18  and  4.20  thru 
4.22.  Direct  comparison  of  the  TSL  procedure  to  the  MEP  is  provided  in 
Figures  4.20  thru  4.22.  Figures  4.18  and  4.20  thru  4.22  are  for  the 
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.19.  Sediment  load  determined  by  MEP  versus 
estimated  sediment  load  for  the 
Niobrara  River. 
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Figure  4.21.  Total  sediment  load  determined  by  MEP  and 
TSL  versus  streamflow  for  section  C  of  the 
Middle  Loup  River. 
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Niobrara  River,  Fivemile  Creek,  the  Middle  Loup  River,  and  the 
Rio  Grande  conveyance  channel,  respectively. 

The  MRP  and  the  observed  sediment  loads  versus  streamflow  are  shown 
in  Figure  4.18  for  the  Niobrara  River.  Comparison  of  Figure  4.18  to 
Figure  4.4  indicates  that  single  MEP  measurements  provide  less  data 
scatter  than  single  TSL  measurements.  However,  comparison  of  Fig¬ 
ure  4.18  to  Figure  4.6  indicates  that  the  use  of  the  TSL  procedure  along 
with  time  integration  will  provide  equally  reliable  sediment  load  mea¬ 
surements  . 

Another  comparison  between  the  MEP  and  the  TSL  procedure  using  the 
Niobrara  River  data  is  made  using  Figures  4.5,  4.6,  and  4.19.  Fig¬ 
ure  4.19  provides  a  comparison  of  the  MEP  sediment  load  measurements  to 
Equation  4.1  (the  regressed  function  selected  to  represent  the  observed 
sediment  loads).  In  Figure  4.19,  89  percent  of  the  MEP  sediment  loads 
are  within  ±50  percent  of  the  loads  determined  by  Equation  4.1.  Fig¬ 
ure  4.5  shows  73  percent  of  the  TSL  sediment  loads  to  be  within  ±50  per¬ 
cent  of  Equation  4.1.  In  Figure  4.6,  100  percent  of  the  averaged  TSL 
sediment  loads  are  within  ±50  percent  of  Equation  4.1.  Comparison  of 
these  statistics  indicates  greater  reliability  for  a  single  MEP  measure¬ 
ment  than  a  single  TSL  measurement.  However,  when  the  average  of 
multiple  TSL  measurements  is  considered,  as  in  Figure  4.6,  greater  re¬ 
liability  is  given  to  the  average  of  multiple  TSL  measurements  than  to  a 
single  MEP  measurement. 

The  MEP,  the  TSL,  and  the  observed  sediment  loads  versus  streamflow 
are  directly  compared  for  Fivemile  Creek  in  Figure  4.20.  Figure  4.20 
indicates  good  agreement  between  the  measurements  in  the  lower  discharge 
range  and  the  majority  of  the  measurements  in  the  higher  discharge 
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range.  However,  in  Che  higher  discharge  range,  selected  MEP  measure¬ 
ments  are  larger  and  selected  TSL  measurements  are  smaller  than  the 
observed  sediment  loads.  The  variation  in  the  TSL  measurements  from  the 
observed  measurements  were  discussed  in  Section  4.3.2.  Burkham  and 
Dawdy  (1980)  suggest  that  the  overestimates  of  sediment  load  by  the  MEP 
may  be  related  to  an  underestimate  of  the  roughness  element  height  kg  . 
However,  no  information  on  the  type  of  bedforms  present  during  sampling 
is  available  to  further  assess  this  explanation  for  the  overestimates. 
Underestimates  of  kg  will  lead  to  an  Improper  calculation  of  bed  load, 
suspended-sediment  concentration  distribution,  and  velocity  profile. 
Further  assessment  between  the  MEP  and  the  TSL  measurements  for  Flvemile 
Creek  is  not  possible  since  concurrent  depth-integrated  and  point 
samples  were  not  recorded. 

The  MEP,  the  TSL,  and  the  observed  sediment  loads  versus  streamflow 
are  directly  compared  for  the  Middle  Loup  River  in  Figure  4.21.  Fig¬ 
ure  4.21  indicates  good  agreement  between  all  of  the  measurements. 
Additional  comparisons  between  the  MEP  and  the  TSL  measurements  for  the 
Middle  Loup  River  were  not  assessed  since  depth-integrated  and  point- 
source  suspended-sediment  samples  were  not  concurrently  collected. 

The  MEP,  the  TSL,  and  the  observed  sediment  loads  versus  streamflow 
are  directly  compared  for  the  Rio  Grande  conveyance  channel  in  Fig¬ 
ure  4.22.  Figure  4.22  Indicates  good  agreement  for  the  majority  of  the 
measurements.  Some  of  the  MEP  sediment  loads  are  lower  than  the 
observed  sediment  loads  for  the  larger  sampled  streamflow.  Assuming 
these  measurements  to  be  correct,  greater  reliability  may  be  placed  upon 
the  lower  TSL  measurements  plotted  in  Figure  4.16. 
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Overall,  comparisons  of  the  MEP,  the  TSL  procedure,  and  the  ob¬ 
served  sediment  load  data  show  similar  results.  The  Niobrara  River  data 
indicate  that  the  MEP  may  provide  greater  reliability  with  a  single 
depth-integrated  suspended-sediment  sample  than  the  TSL  procedure  with  a 
single  point  suspended-sediment  sample.  However,  the  Fivemile  Creek, 
the  Middle  Loup  River,  and  the  Rio  Grande  conveyance  channel  data 
indicate  a  single  TSL  measurement  to  be  as  reliable  as  a  single  MEP  mea¬ 
surement.  The  reliability  of  the  TSL  procedure  is  improved  when 
multiple-point  samples  can  be  collected  and  Integrated  into  the  analy¬ 
sis.  The  similarity  of  the  TSL,  the  MEP,  and  the  observed  sediment 
loads  shows  promise  for  use  of  the  TSL  procedure  to  measure  total  sedi¬ 
ment  load. 

4.5  APPLICATION  TO  THE  NEWLY 

ACQUIRED  GOODWIN  CREEK  DATA 

The  final  objective  of  this  study  was  to  apply  the  TSL  procedure  to 
newly  acquired  data.  The  new  data  was  provided  by  the  US  Department  of 
Agriculture-Agricultural  Research  Service  (USDA-ARS) ,  Sedimentation 
Laboratory,  Oxford,  Mississippi,  and  was  collected  at  the  Goodwin  Creek, 
Mississippi,  research  catchment  (Seely,  Grissinger,  and  Little  1981, 

US  Army  Corps  of  Engineers  1981,  and  Bowie  and  Sansom  1986).  The 
Goodwin  Creek  research  catchment  has  14  gaging  stations  within  an 
8.2-square-mile  drainage  area.  The  gaging  stations  have  supercritical- 
flow  flumes  with  instrumentation  to  record  water  level,  water  tempera¬ 
ture,  and  precipitation.  Each  gaging  station  has  an  automatic  pumping- 
type  suspended-sediment  sampler.  The  outlet  of  the  flumes  is  designed 
with  a  flow  overfall  which  permits  total  sediment  load  measurements  with 
the  absence  of  an  unmeasured  zone  in  the  measured  vertical.  The  total 
load  measurements  are  made  with  a  depth-integrated  sampler. 


The  USDA-ARS  Sedimentation  Laboratory  provided  detailed  data  from 
Goodwin  Creek,  Station  2,  for  the  assessment  of  the  TSL  procedure  for  a 
storm  event  occurring  on  August  27,  1986.  Unlike  other  sediment  data, 
the  August  27,  1986,  data  set  for  Goodwin  Creek  provides  concurrent  and 
consecutive  point-source  suspended-sediment  samples  and  total  sediment 
load  measurements  for  a  streamflow  event.  The  data  from  the  streamflov 
event  that  was  used  with  the  TSL  procedure  are  provided  in  Appendix  D. 
The  assessment  of  TSL  with  these  data  is  accomplished  using  Fig¬ 
ures  4.23-4.27.  Total  sand  load  is  assessed  in  these  figures  rather 
than  total  sediment  load,  because  Goodwin  Creek's  bed  material  is  ap¬ 
proximately  80  percent  gravel,  which  is  not  measured  by  the  depth- 
integrated  and  automatic  pumping-type  suspended  sediment  samplers. 

Figure  4.23  shows  the  hydrograph  for  the  AuguBt  27,  1986,  stream- 
flow  event.  Figure  4.24  presents  both  the  total  sand  load  observed  at 
the  flume's  outlet  using  a  depth- integrated  sampler  and  the  total  sand 
load  determined  by  the  TSL  procedure  versus  the  streamflow.  Figure  4.24 
does  not  indicate  a  smooth  function  to  describe  the  relation  between 
sediment  transport  and  streamflow;  therefore,  data  plots  which  show  the 
sediment  load  through  time  for  the  streamflow  event  are  presented  in 
Figures  4.25,  4.26,  and  4.27.  In  these  figures,  greater  fluctuations 
exist  between  consecutively  observed  sand  loads  than  between  consecutive 
TSL  procedure  sand  loads.  Figure  4.27  shows  the  magnitude  of  the 
fluctuations  for  the  TSL  procedure  sand  loads  to  be  similar  to  or  to 
occur  within  the  fluctuations  of  the  observed  sand  loads.  This  obser¬ 
vation  may  be  attributed  to  differences  in  the  two  sampling  techniques. 
Samples  for  the  observed  total  sand  loads  are  collected  with  a  depth- 
integrated  sampler  being  lowered  and  raised  through  the  entire  vertical 
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Figure  4.24.  Total  sand  load  versus  streamflov  for  the 

August  27,  1986,  streanflow  event  at  Goodwin 
Creek,  Mississippi,  Station  2. 
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Figure  4.26.  Total  sand  load  measured  by  the  TSL  procedure 
versus  streamflow  for  the  August  27,  1986, 
streamflow  event  at  Goodwin  Creek,  Mississippi, 
Station  2. 
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Figure  4.27.  Total  sand  load  versus  time  for  the  August  27, 
1986,  streamflov  event  at  Goodwin  Creek, 
Mississippi,  Station  2.  Suspended-sediment 
point  samples  are  from  20  to  45  percent  of 
the  flow  depth  above  the  channel  bed. 
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flow  field.  The  TSL  suspended-sediment  samples  are  collected  at  one 
point. 

4.6  DISCUSSION 

The  TSL  procedure  was  compared  to  total  sediment  load  field  data 
and  the  MEP.  The  assessed  data  were  collected  in  small  streams  with 
flow  depths  of  less  than  5  ft.  Four  of  the  assessed  data  sets  were  for 
the  Niobrara  River  at  Cody,  Nebraska;  Fivemile  Creek  at  Shoshoni, 
Wyoming;  the  Middle  Loup  River  near  Dunning,  Nebraska;  and  the  Rio 
Grande  conveyance  channel  near  Bernardo,  New  Mexico.  These  data  sets 
had  previously  been  used  in  the  development  and  evaluation  of  the  MEP. 

A  newly  acquired  data  set  collected  at  Goodwin  Creek,  Mississippi,  was 
also  assessed. 

Comparison  of  the  TSL  procedure  to  field  data  and  the  MEP  highlight 
the  use  of  the  TSL  procedure  as  an  alternative  to  depth-integrated 
suspended-sediment  sampling.  Statistics  for  the  Niobrara  River  data 
indicate  that  the  MEP  procedure  may  provide  greater  reliability  with  a 
single  depth- integrated  suspended-sediment  sample  than  the  TSL  procedure 
with  a  single  point  suspended-sediment  sample.  The  Niobrara  River  data 
showed  89  percent  of  the  the  MEP  measurements  and  76  percent  of  the  TSL 
procedure  measurements  to  be  less  than  a  50  percent  deviation  from  sedi¬ 
ment  load  estimated  by  Equation  4.1.  Further  assessment  of  the  Niobrara 
River  data  Indicates  that  the  average  of  multiple  TSL  measurements,  as 
in  Figure  4.6,  is  as  reliable  as  a  single  MEP  measurement. 

The  Fivemile  Creek,  the  Middle  Loup  River,  and  the  Rio  Grande  con¬ 
veyance  channel  data  Indicate  a  single  TSL  measurement  is  as  reliable  as 
a  single  MEP  measurement.  Comparisons  on  the  percent  differences 
between  the  MEP  and  the  TSL  procedure  for  these  data  were  not  assessed 
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since  concurrent  measurements  of  depth-integrated  and  point  suspended- 
sediment  samples  were  not  reported. 

The  tested  data  show  that  the  measurement  procedure  is  improved 
when  point  suspended-sediment  samples  are  obtained  in  the  lower  flow 
depth  but  not  too  close  to  the  channel  bed.  Lane  and  Kalinske  (1941) 
recommended  that  point  samples  should  be  collected  in  the  lower  half¬ 
depth.  The  TSL  procedure  was  tested  and  found  to  be  reliable  using 
point  samples  collected  in  the  lower  20  to  45  percent  flow  depth.  Point 
samples  collected  at  locations  further  from  the  channel  bed  may  provide 
less  reliable  measurements  due  to  deviations  in  the  theoretical  sediment 
concentration  and  flow  velocity  profiles.  These  deviations  may  result 
from  macroturbulence  and  are  not  numerically  significant  to  the  TSL  pro¬ 
cedure  when  the  point  samples  are  collected  in  the  lower  flow  depth. 
Point  samples  collected  close  to  the  channel  bed  may  exhibit  greater 
variance  than  other  point  samples  due  to  the  same  processes  that  influ¬ 
ence  the  stochastic  nature  of  bed  load  transport.  In  particular,  point 
samples  may  not  be  representative  if  they  are  collected  near  bed  forms 
of  a  significant  height. 

The  Goodwin  Creek  data  provided  consecutive  and  concurrent  depth- 
integrated  and  point  suspended-sediment  samples.  These  data  demonstrate 
that  both  depth-integrated  measurements  and  TSL  procedure  measurements 
can  yield  scatter  in  sediment  load  through  time.  Comparison  of  the 
total  sand  load  observed  with  the  depth- integrated  sampler  to  the  total 
sand  load  measured  using  the  TSL  procedure  demonstrates  that  the  TSL 
procedure  provides  a  reliable  measurement  of  sediment  load. 

The  assessments  of  the  reviewed  data  sets  highlight  advantages  for 
selecting  the  TSL  procedure  over  the  MEP.  Advantages  of  the  TSL 
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procedure  Include  realistic  sediment  load  measurements,  greater  control 
In  obtaining  suspended-sediment  samples,  automatic  sampling,  the  collec¬ 
tion  of  near-continuous  measurements,  and  time  integration  for  determin¬ 
ing  the  true  sediment  transport.  Also,  the  use  of  point-source 
suspended-sediment  samples  requires  less  time  in  the  field  and  lower 
sediment  sampling  costs. 


CHAPTER  5 


SUMMARY,  CONCLUSION,  AND  RECOMMENDATIONS 

5.1  SUMMARY 

The  main  objective  of  this  study  was  to  develop  a  procedure  for 
total  sediment  load  measurement  using  a  point-source  suspended-sediment 
sample.  The  procedure  was  developed  and  compared  to  field  data  and  the 
modified  Einstein  procedure  (MEP) .  The  computational  step  solution  and 
computer  code  listing  for  the  proposed  procedure  are  provided  in 
Appendices  A  and  B,  respectively.  The  proposed  procedure  and  computer 
code  are  identified  as  TSL  (Total  Sediment  Load) . 

The  TSL  procedure  applies  Einstein's  (1950)  modification  to 
Keulegan's  (1938)  mean  and  point  velocity  relations,  the  suspended- 
sediment  concentration  distribution  equation  developed  by  Rouse  (1937), 
and  Burkham  and  Dawdy's  (1980)  modified  form  of  Einstein's  bed  load 
function.  The  data  required  for  application  of  the  TSL  procedure 
include  descriptors  of  the  streamflow,  the  channel  geometry,  the  bed 
material,  and  the  point-source  suspended-sediment  sample. 

The  TSL  procedure  was  evaluated  using  field  data  from  the  Niobrara 
River  at  Cody,  Nebraska;  Flvemile  Creek  at  Shoshonl,  Wyoming;  the  Middle 
Loup  River  near  Dunning,  Nebraska;  the  Rio  Grande  conveyance  channel 
near  Bernardo,  New  Mexico;  and  Goodwin  Creek  near  Oxford,  Mississippi. 
The  field  data  enabled  comparisons  of  the  TSL  procedure,  the  MEP,  and 
field  measurements  of  total  sediment  load.  The  similarity  in  the 
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results  of  the  comparisons  highlights  the  use  of  the  TSL  procedure  as  an 
alternative  to  depth-integrated  suspended-sediment  sampling. 

5.2  CONCLUSIONS 

The  TSL  procedure  provides  total  sediment  load  measurements  that 
are  similar  in  magnitude  to  both  field  and  MEP  measurements.  The 
Niobrara  River  data  indicate  that  the  MEP  procedure  may  provide  greater 
reliability  with  a  single  depth-integrated  suspended-sediment  sample 
than  the  TSL  procedure  with  a  single  point  suspended-sediment  sample. 
However,  further  assessment  of  the  Niobrara  River  data  indicates  that 
the  average  of  multiple  TSL  measurements  is  as  reliable  as  a  single  MEP 
measurement.  The  Fivemile  Creek,  the  Middle  Loup  River,  and  the 
Rio  Grande  conveyance  channel  data  indicate  a  single  TSL  measurement  is 
as  reliable  as  a  single  MEP  measurement. 

The  tested  data  show  that  the  TSL  measurement  procedure  is  Improved 
when  point  suspended-sediment  samples  are  obtained  in  the  lower  flow 
depth.  Lane  and  Kalinske  (1941)  recommended  that  point  samples  should 
be  collected  in  the  lower  half-depth.  The  TSL  procedure  was  tested  and 
found  to  be  reliable  using  point  samples  collected  in  the  lower  20  to 
45  percent  flow  depth. 

The  Goodwin  Creek  data  provided  consecutive  and  concurrent  depth- 
integrated  and  point  suspended-sediment  samples.  These  data  demonstrate 
that  both  depth-integrated  measurements  and  TSL  procedure  measurements 
can  yield  scatter  in  sediment  load  through  time.  Comparison  of  the 
total  sand  load  observed  with  the  depth-integrated  sampler  to  the  total 
sand  load  measured  using  the  TSL  procedure  demonstrates  that  the  TSL 
procedure  provides  a  reliable  measurement  of  sediment  load. 
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The  assessments  of  the  reviewed  data  sets  highlight  advantages  for 
selecting  the  TSL  procedure  over  the  MEP.  Advantages  of  the  TSL  proce¬ 
dure  include  realistic  sediment  load  measurements,  greater  control  in 
obtaining  suspended-sediment  samples,  automatic  sampling,  the  collection 
of  near-continuous  measurements,  and  time  integration  for  determining 
the  true  sediment  transport.  Also,  the  use  of  point-source  suspended- 
sediment  samples  requires  less  time  in  the  field  and  lower  sediment 
sampling  costs. 

5 . 3  RECOMMENDATIONS  FOR 

FURTHER  STUDIES 

The  data  used  to  compare  the  TSL  procedure  to  turbulence  flume 
total  sediment  load  measurements  and  the  MEP  were  not  from  concurrently 
collected  data  sets.  The  exception  was  concurrently  collected  data  at 
the  Goodwin  Creek  research  catchment.  Additional  data  sets  which  pro¬ 
vide  concurrent  sampling  of  point-source  suspended-sediment  samples 
should  be  collected  and  assessed.  The  additional  data  sets  will  assist 
in  learning  more  about  field  procedures  for  sediment  measurement  and  the 
sediment  transport  processes  that  each  sediment  sampling  procedure  is 
measuring.  For  example,  a  depth-integrated  sample  measures  a  flow 
column  over  time  while  a  point  sample  can  measure  a  point  over  time. 
Further  studies  may  confirm  that  point  sampling  with  application  of  the 
TSL  procedure  and  time  integration  will  provide  more  realistic  measure¬ 
ments  of  sediment  transport. 

The  TSL  procedure  now  permits  researchers  to  consider  continuous 
sampling.  The  research  community  is  now  encouraged  to  develop  a  point 
sampler  which  can  continuously  collect  or  record  data  representing 
suspended-sediment  size  analyses  and  concentrations  at  a  point.  The 
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proposed  sampling  device  along  with  the  TSL  procedure  will  permit  con¬ 
tinuous  measurement  of  total  sediment  load  through  time. 

Further  assessment  should  be  made  on  the  location  of  the  sediment 
sampler  in  the  vertical  and  in  the  cross  section.  The  assessment  will 
require  continuous  total  sediment  load  sampling  and  a  near-constant 
sediment  load  and  streamflow.  The  present  recommendation  is  to  collect 
a  suspended-sediment  sample  in  the  lower  flow  depths  (the  lower  20  to 
45  percent  flow  depth  was  assessed  in  this  study)  in  a  vertical  section 
representing  the  cross  section.  Additional  studies  may  consider  multi¬ 
ple  point  samples  across  a  cross  section  with  one  sampler  per  vertical 
section. 

Studies  designed  to  assess  and  improve  the  determination  of  the 
vertical  velocity  distribution,  the  sediment  concentration  distribution, 
and  the  bed  load  discharge  will  be  useful  for  further  modifications  to 
the  TSL  procedure  and  the  general  understanding  of  sediment  transport. 
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APPENDIX  A 


STEP  SOLUTION 


The  following  computational  procedure  may  be  used  to  determine  the 

total  sediment  load  using  a  point-source  suspended-sediment  sample. 

Step  1.  Identify  the  basic  data:  Q  ,  streamflow  for  the  sampled  cross 

section;  W  ,  width  of  the  cross  section;  A£  ,  cross- 

sectional  area;  R  ,  hydraulic  radius;  d  ,  flow  depth  of  the 

s 

sampled  vertical;  y^  ,  sample  height  above  the  channel  bed; 

C  ,  concentration  of  the  suspended-sediment  sample,  mg/liter; 
P 

i  ,  fraction  of  bed  material  in  each  size  range;  i  ,  frac- 

D  S 

tion  of  suspended  material  in  each  size  range;  , 

and  Dg^  ,  particle  sizes  at  which  35,  65,  and  84  percent  by 
weight,  respectively,  of  the  bed  material  is  finer,  mm  ;  and 
v  ,  the  kinematic  viscosity  of  the  water.  (Except  when  noted 
otherwise,  the  English  system  of  units  is  used  throughout  the 
computation. ) 

Step  2.  Compute  the  mean  flow  depth  d  and  the  average  flow  velocity 
V  using 


d 


A.  1 


and 


V  - 


A. 2 
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Step  3. 


Step  4. 


Step  5. 


Step  6. 


90 

Determine  the  effective  grain  roughness  height  kg  .  The 
computer  program  TSL  (Total  Sediment  Load)  provides  a  default 
of  kg  ■  3.3D84  as  proposed  by  Limerinos  (1970).  Otherwise, 
the  program  user  may  supply  a  Dg4  multiplier  for  kg  . 
Compute  the  shear  velocity  due  to  grain  roughness  using 

the  modified  form  of  Equation  3.3, 


A. 3 


where  x'  is  a  function  of  kg/6'  which  is  presented  in 
Figure  4  of  Einstein  (1950)  and  6'  -  (11.6v)/U^  . 

Compute  the  RS’  product  from 


RS' 


A. 4 


where  S'  is  the  energy  grade  line  slope  due  to  grain 
roughness  and  g  is  the  acceleration  due  to  gravity. 

Determine  the  effective  roughness  height  kg  due  to  grair  and 
form  roughness.  The  computer  program  TSL  provides  a  default  of 
kg  -  3.3Dg4  as  proposed  by  Limerinos  (1970).  Other  available 
user  options  Include  the  procedures  of  Einstein  (1950), 

Brownlie  (1981,  1983),  or  a  user-supplied  Dg4  multiplier. 

When  either  the  Brownlie  or  the  Einstein  procedure  is  applied. 


the  friction  factor  f  is  determined  and  the  relation 


where  is  the  shear  velocity,  is  applied  with  Equation  2.19 


Step  7. 


Step  8. 


for  the  determination  of  k 

s 

Compute  the  shear  velocity  due  to  grain  and  form  roughness 
using  Equation  2.19  modified  as 


A.  6 


where  x  is  a  function  of  kg/6  which  is  presented  in 
Figure  4  of  Einstein  (1950)  and  6  ■  11.6v/UA  . 

Compute  the  intensity  of  shear  for  each  bed  material  size 

interval ,  where 


l.escsY)^ 


A. 7 


in  which 


8^  ■  log  10.6 


8  .  -  log  (10. 6x) 


6  +  C  D 

(£Y)ow  *  - i — ~  when  D  <  &  +  C  D  and 

KM  K  S 

S 

*  t—  when  D  i  6  +  C  D 
k  s 

s 

and 

C  -  -0.62  +  3.12(D,t.D,,.)1/2 

S  03  J3 

(The  representative  grain  size  D  for  each  sediment  size 
interval  is  the  geometric  mean  of  the  interval  limits.) 

Step  9.  Determine  the  transport  intensity  for  each  bed  material 

size  interval  using  Figure  10  from  Einstein  (1950)  and  i|^ 
from  step  8. 

Step  10.  Compute  the  bed  load  discharge  i_q_  for  a  unit  width  of 

D  D 

flow  for  each  bed  material  size  interval  using  Equation  2.28. 

2 

With  g  -  32.2  ft/sec  ,  the  density  of  solids  pg 
3 

*  5.17  slugs/ft  ,  and  the  specific  gravity  of  the  sediment 
Sg  ■  2.65  ,  Equation  2.28  becomes 

1^3  -  1200ibD3/\  A. 8 

where  is  the  fraction  of  bed  material  in  a  given  size 

range . 

Step  11.  Compute  the  bed  load  concentration,  C.  ,  for  each  bed 

o 

material  size  fraction  using 
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„  ,  _J_  iBqB 
Jb  11.6  2DUJ 


A. 9 


where  11. 6U^  is  equlvalenc  to  the  race  of  sediment  transport 
in  the  bed  layer  as  proposed  by  Einstein  (1950). 

Step  12.  Compare  the  concentration  of  the  suspended  sediment  sample 
to  the  bed  load  concentration  .  If  »  set 

-  Cp  and  recompute  i^q^  •  This  limits  the  suspended- 
sediment  concentration  at  less  than  or  equal  to  the  bed  load 
concentration. 

Step  13.  Determine  the  lower  limit  of  the  suspended  load  zone.  This 

limit  is  set  at  the  largest  of  two  possible  heights  above  the 
channel  bed.  The  first  limiting  height  is  y^  *  2D 
(Einstein’s  suggestion  for  the  height  of  the  bed  load  zone) 
and  the  second  limiting  height  is  where  »  0  according  to 
Equation  2.18,  where  is  the  average  flow  velocity  at  a 

distance  y  from  the  channel  bed  or 


30. 2x 


A.  10 


Step  14.  Determine  z  for  each  measured  suspended  sediment  size 
interval  using  Equation  2.26  from  which 
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Step  15.  Develop  a  linear  regressed  relationship  for  z 


z  ■  djW 


A. 12 


where  and  a^  are  regression  coefficients  and  w  Is  the 

particle  fall  velocity  for  the  geometric  mean  of  each  sampled 
sediment  size  interval  as  determined  using  Ruby’s  equation 
(Ruby  1933) 


' 

»  * 

2  t  36v2 

1/2 

2 

36v 

1/2 

3  g(Ss  -  1)D3 

g(S8  -  1)D3 

b 

- 

•  » 

4 

A.  13 


[g(Sa  -  1)D] 


1/2 


Step  16.  Recompute  C.  for  each  sampled  suspended  sediment  size 

D 

fraction  with  Equation  2.26  modified  as 


d  -  y„  y. 


A.  14 


£  __t 


yv  d  - 


wherein  y.  is  set  for  each  sediment  size  Interval  according 

D 

to  step  13  and  z  is  determined  using  Equation  A. 12. 


Step  17.  Recomputed  i_q_  for  each  sampled  suspended  sediment  size 

D  D 


fraction,  where 


Vb  -  Sn-6(2D)Di 


A. 15 


The  height  of  the  bed  load  zone  is  set  at  2D  for  all 
sediment  size  fractions  for  the  numerical  calculations  to 
agree  with  Einstein's  (1950)  formulation  of  bed  load 
transport. 

Step  18.  Compute  the  suspended  sediment  transport  i^q^  for  each 
sampled  size  fraction  using 

iaq8  "  +  *2^  A*16 

with  P  given  by  Equation  2.52  Ij  given  by  Equation  2.48, 
and  1 2  given  by  Equation  2.49. 

Step  19.  Compute  the  total  sediment  load  through  a  cross  section 
using 

1TqT  “  iBqB  +  1sqs  A,1? 

and 

Qt  -  43.2wEiTqT  A. 18 

where  43. 2W  is  the  stream  width  multiplied  by  a  constant  for 
converting  pounds  per  second  to  tons  per  day. 


APPENDIX  B 


COMPUTER  PROGRAM  LISTING 


The  following  computer  program  was  written  for  use  on  a  personal 
computer.  The  program  Is  written  In  Fortran  77  and  can  be  run  In 
executlble  form  on  an  IBM  compatible  personal  computer  with  128K  and  a 
single  floppy  disk  drive.  During  program  execution,  the  user  Is  prompted 
for  Input  and  output  data  file  names. 


PROGRAM  TSL 
C  VERSION;  9-29-87 
CHARACTER* 8  FL 

REAL  IBQB(15) ,ISQS(15) ,ITQT(15) ,KS,KSP 
DIMENSION  D(15) ,FB(15) ,FS(15) ,PB(15) ,PS(15) ,CA(15) , 
1CB(15) ,PSI(15) , PHI (15) ,GM(15) ,Z(15) ,Z1(15) ,W(15) , 
2A(15) 

C 

C 

C********************************************************** 

C  INPUT  VARIABLES 

C********************************************************** 

C 

C  AREA  -  FLOW  AREA  OF  MEASURED  CROSS  SECTION,  FT** 2 
C  CON  -  CONCENTRATION  OF  SUSPENDED  SEDIMENT  SAMPLE, 

C  MG/LITER 

C  CS  -  COEFFICIENT  USED  TO  COMPUTE  BED  LOAD  WHICH  IS  A 
C  FUNCTION  OF  MEDIAN  DIAMETER 

C  IF  CS-O,  CS  WILL  BE  COMPUTED 

CD-  SEDIMENT  SIZE  FRACTION  DIVISIONS,  MM 
C  DKS  -  EFFECTIVE  FLOW  ROUGHNESS  HEIGHT  MULTIPLIER 
C  (KS-DKS*D84)  IF  DRS-0,  KS  WILL  BE  COMPUTED  (USE 

C  OF  THIS  PARAMETER  OVERRIDES  THE  USE  OF  IKS) 

C  DKSP  -  EFFECTIVE  GRAIN  ROUGHNESS  HEIGHT  MULTIPLIER 

C  (KSP-DKSP*D84) 

C  IF  DKSP-0,  DKSP  DEFAULTS  TO  3.3 

C  DSV  -  DEPTH  OF  SAMPLING  VERTICAL,  FT 
C  IKS  -  IDENTIFIER  FOR  FRICTION  FACTOR  COMPUTATION 
C  METHOD 


C  IKS-0  FOR  LIMERINOS 

C  IKS-1  FOR  BROWNLIE 

C  IKS-2  FOR  EINSTEIN 

C  ISTA  -  STATION  OR  SAMPLE  NUMBER 

C  ND  NUMBER  OF  SIZE  FRACTIONS  (MAXIMUM  -  12) 

C  NST  -  NUMBER  OF  TOTAL  LOAD  SAMPLES 

C  PB  PERCENT  OF  BED  MATERIAL  FINER  THAN  D 

C  PS  PERCENT  OF  SUSPENDED  MATERIAL  FINER  THAN  D 

C  Q  STREAMFLOW,  CFS 

C  R  HYDRAULIC  RADIUS  OF  MEASURED  CROSS  SECTION,  FT 
C  T  WATER  TEMPERATURE,  DEGREES  FAHRENHEIT 

C  WD  CROSS  SECTION  WIDTH,  FT 

C  YS  SAMPLE  HEIGHT  ABOVE  CHANNEL  BED,  FT 
C 
C 

C********************************************************** 

C*  INPUT  DATA  FILE  FORMAT  REQUIREMENTS 

C******************************** ******* ******************* 

c 

C  READ(11 ,230)  ND 

C 

C  DO  5  1-1  ,ND+1 

C  5  READ(11 ,210)  D(I) 

C 

C  READ(11 ,230)  NST, IKS 

C 

C  DO  900  I ST- 1, NST 

C  READ(11 ,240)  Q,WD, AREA, R.DSV.YS.T, ISTA 

C  READ(1 1,255)  (PB(J) , J-l ,ND+1) 

C  READ(11,255)  (PS(J)  ,  J-l  ,NDH) 

C  READ(11 ,260)  CON,CS,DKS,DKSP 

C  900  CONTINUE 
C 

C  210  F0RMAT(F10.4) 

C  230  F0RMAT(2I5) 

C  240  FORMAT (7F1 0.2,15) 

C  255  FORMAT (15F5.0) 

C  260  FORMAT(4F10.0) 

C 

C********************************************************** 

c 

WRITE (*,*)’ ENTER  INPUT  DATA  FILE  NAME’ 

READ (*,200)  FL 
200  FORMAT(A8) 

OPEN ( 1 1 , FILE-FL ) 

WRITE (*,*)' ENTER  OUTPUT  FILE  NAME' 

READ (*,200)  FL 

OPEN (12, FILE-FL , STATUS- ' NEW' ) 

WRITE(*,*)  ’ENTER  GRAPHICS  OUTPUT  FILE  NAME’ 
READ(*,200)  FL 

OPEN (13, FILE-FL , STATUS- ’ NEW ' ) 


c 

c********************************************************** 

C*  READ  INPUT  DATA  FILE  /  SIZE  FRACTION  DIVISIONS 

C********************************************************** 

C 

READ(1 1,230)  ND 
DO  5  I-l.ND+l 
READ(11 ,210)  D(I) 

210  FORMAT (F10. 4) 

5  CONTINUE 
C 

C********************************************************** 

C*  DETERMINE  EACH  SIZE  RANGE’S  GEOMETRIC  MEAN,  FT 

C********************************************************** 

C 

DO  10  I-l.ND 

10  GM(I)*(D(I)*D(I+1) )**0. 5/304 .8 
C 

C********************************************************** 

C*  READ  AND  ANALYZE  EACH  DATA  SET 

c 

READ(11,230)  NST , IKS 
230  F0RMAT(2I5) 

WRITE(13,205)IKS 

205  FORMAT ( ’ IKS  -  ’ ,I5,/,11X,'Q' ,7X,'TL' ,8X,'TF' ,8X,'TS’ , 
19X,’TG’  , 7X , ’ TFS ' ,3X,'YS/DSV' ,1X,’DSV/D’ ,1X,'STA' ,3X, 
2'D50G’ ,7X,’CS’ ,/) 

DO  900  IST-l.NST 
DO  15  J-l.ND 
Z(J)-0.0 
15  Z1(J)«0.0 

WRITE(*,235)  1ST, NST 

235  FORMAT(////,’  DATA  SET  ',13,'  OF  ’,13,//) 

READ(11 ,240)  Q,WD,AREA,R,DSV,YS ,T,ISTA 
240  FORMAT(7F10. 2 ,15) 

DEPTH-AREA/WD 
V-Q/AREA 
Xl-YS/DSV 
X2-DSV /DEPTH 

WRITE(12 ,250)  Q, V,VTO, DEPTH, DSV, YS ,T, IKS 
250  FORMAT(////,’TSL  VERSION  9-29-87',///, 

1 ' STREAMFLOW , CFS  -  ' ,£10.3,/, 'MEAN  VELOCITY,  FPS  -  ' , 
2E10. 3, /.'CHANNEL  WIDTH,  FT  -  '.E10.3,/, 

3 'MEAN  FLOW  DEPTH,  FT  -  'E10.3,/, 

4 'DEPTH  OF  SAMPLING  VERTICAL,  FT  -  ’.E10.3,/, 

5' SAMPLE  HEIGHT  ABOVE  CHANNEL  BED,  FT  -  '.E10.3,/, 

6 'WATER  TEMPERATURE,  DEGREES  FAHRENHEIT  -  '.E10.3,/, 

7 'IKS  -  ' ,15) 

READ(11 ,255)  (PB(J) , J-l ,ND+1) 

READ (11 ,255)  (PS(J) ,J-1 ,ND+1) 

READ(11 ,260)  CON , CS , DKS , DKSP 
255  FORMAT (15F5.0) 

260  FORMAT (4F10.0) 


c********************************************************** 
C  DETERMINE  D16.D35 ,D50,D65,  AND  D84;  MM 

C********************************************************** 

C 

DO  20  J-l.ND 

FB(J)-(PB(J+1)-PB(J))/100.0 
FS(J)«(PS(J+1)-PS(J))/100.0 
IF(FB(J) .LE.O.O)  FB(J)-0.0 
IF(FS(J) .LE.O.O)  FS(J)-0.0 
20  CONTINUE 
P16L-0.0 
DO  25  J-I.ND+I 

IF(P16L.GT.O.O.AND.PB(J).LE.O.O)GO  TO  25 

IF(PB(J) .LE.16.0)  THEN 

PI6L-PB(J) 

D16L-D(J) 

END  IF 

IF(PB(J) .LE.35.0)  THEN 
P35L-PB(J) 

D35L-D(J) 

ENDIF 

IF(PB(J) .LE.50.0)  THEN 
P50L-PB(J) 

D50L-D(J) 

ENDIF 

IF(PB(J) .LE.65.0)  THEN 
P65L-PB(J) 

D65L-D(J) 

ENDIF 

IF(PB(J) .LE.84.0)  THEN 
P84L«PB(J) 

D84L-D(J) 

ENDIF 

25  CONTINUE 

DO  30  J-ND+I ,1,-1 
IF(PB(J) .GT.16.0)  THEN 
P16H-PB(J) 

D16H-D(J) 

ENDIF 

IF(PB(J) .GT.35.0)  THEN 
P35H-PB(J) 

D35H-DCJ) 

ENDIF 

IF(PB(J) .GT.50.0)  THEN 
P50H-PB(J) 

D50H-D(J) 

ENDIF 

IF(PB(J) .GT.65.0)  THEN 
P65H-PB(J) 

D65H-D(J) 

ENDIF 

IF(PB(J) .GT.84.0)  THEN 
P84H-PB(J) 


D84H-D(J) 

ENDIF 

30  CONTINUE 
AB-P16H-P16L 

B-ALOG 1 0 (D 1 6H) -ALOG 1 0 (D 1 6L) 

SM-AB/B 

AB- ( 16 . 0-P 1 6L) / SM+ALOG 1 0 (D 1 6L) 

D16-10.0**AB 

AB-P35H-P35L 

B»ALOG10(D35H)-ALOG10(D35L) 

SM-AB/B 

AB-(35.0-P35L)/SM+ALOG10(D35L) 

D35-10.0**AB 

AB-P50H-P50L 

B-ALOG10(D50H)-ALOG10(D50L) 

SM-AB/B 

AB-(50.0-P50L) / SM+ALOG 10(D50L) 

D50-10 .0**AB 
AB-P65H-P65L 

B- ALOG 1 0 (D6  5H ) - ALOG 1 0 (D  6  5L ) 

SM-AB/B 

AB-(65.0-P65L) / SM+ALOG 10(D65L) 

D65-10.0**AB 

AB-P84H-P84L 

B-ALOG10(D84H)-ALOG10(D84L) 

SM-AB/B 

AB-(84.0-P84L) / SM+ALOG 10(D 8 4L) 

D84-10.0**AB 
WBITE(12,*) *D16  -  ' ,D16 
WRITE(12,*)’D35  -  ’ ,D35 
WRITE(12,*) 'D50  -  ’ ,D50 
WRITE (12,*)' D65  -  ' ,D65 
WRITE ( 1 2 , *) ' D84  -  *,084 
C 

C*********** *********************************************** 

C*  DETERMINE  THE  SHEAR  VELOCITY  POR  GRAIN  ROUGHNESS, 

C*  FT/ SEC 

C********************************************************** 

c 

VIS-VISCOS(T) 

IF(DKSP) 85 ,85,90 
85  DKSP-3.3 
90  KSP-DKSP*D84 / 304 . 8 
RKS-KSP 

CALL  SHEAR (V.R.RKS, VIS, XP.USP) 

SP-USP**2.0/(32.2*R) 

DKSP»KSP*304.8/D84 
DKSPA-KSP*  304.8 /D6  5 

WRITE(12,*)'KSP  -  ' ,KSP, ’  KSP/D84  -  ' ,DKSP, 

1'  KSP/D65  -  ' ,DKSPA 
WRITE (12,*) *  SP  -  ' ,SP,'  XP  -  ' ,XP 
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C 

c********************************************************** 
C*  DETERMINE  KS,  FT,  EFFECTIVE  FLOW  ROUGHNESS  HEIGHT 
C********************************************************** 

c 

1*0 

FG-V/SQRT ( 1 .65*32 . 2*D50/304 . 8) 

IF(DKS)35,35,33 

C 

C********************************************************** 

C  DETERMINE  FRICTION  FACTOR  USING  THE  ENTERED  DKS 

C********************************************************** 

C 

33  RKS*DKS*D84 / 304 . 8 
KS-RKS 

CALL  SHEAR(V,R,RKS,VIS,X,US) 

F*8 . 0* (US/V) **2 . 0 
GO  TO  80 

35  GO  TO(36, 37 ,42)IKS+1 
C 

C***************** ******* ********************************** 

C  DETERMINE  FRICTION  FACTOR  USING  LIMERINOS  (1970) 

C********************************************************** 

c 

36  DKS-3.3 
GO  TO  33 

C 

C********************************************************** 
C  DETERMINE  FRICTION  FACTOR  USING  BROWNLIE  (1981) 

£***** ****** ************************** ********************* 

C 

37  D50-D50/304.8 
SIG-SQRT(D04/D16) 

S-0.038805*FG**2.572*SIG**0.413/(R/D50)**1.361 
IF(S.LE. 0.006) GO  TO  39 

38  S-0.021425*FG**2. 172*SIG**0. 2785/ (R/D50)**i .304 
GO  TO  41 

39  FGP-1. 74/S**. 3333 
IF(FGP/FG~l .0) 41,41,40 

40  F-0.188925*FG**0.572*SIG**0.413/(R/D50)**0.361 
GO  TO  65 

41  I-I+l 

IF(I.LE. l)GO  TO  38 

F-0 . 104346*FG**0 . 1 72*SIG**0 . 279 / (R/D50) **0. 304 
GO  TO  65 
C 

C********************************************************** 

C  DETERMINE  FRICTION  FACTOR  USING  EINSTEIN  (1950) 

Q************************** ******************************** 

C 

42  RSP«USP**2 .0/32.2 
CHI-1 . 65*D35/ (304 .8*RSP) 

CALL  FIG5(CHI ,USPP) 


USPP-V/USPP 

US-SQRT(USP**2 . 0+USPP**2 . 0) 

F-8.0/(V/US)**2.0 

C 

65  DEL-12.27*R/10.0**(SQRT(8.0/F)/5.75) 

US-SQRT(F/8.0)*V 

DELTA- 1 1 . 6*VIS/US 

X-1.0 

KS-DEL*X 

SIN-KS /DELTA 

XP-FX(SIN) 

IF(XP-X) 70,80,70 
70  KS-20.0*DELTA 
75  SIN-KS /DELTA 
X-FX(SIN) 

KSP-X*DEL 

TEST-ABS (KSP-KS ) /KS 
IF(TEST.LE. 0.001) GO  TO  80 
KS-KSP 
GO  TO  75 

80  DKS-KS/D84*304 . 8 
DKSA-KS/D65*304.8 
S-US**2. 0/(32. 2*R) 

FGP-1. 74/S**. 3333 

WRITE (12,*) 'KS  -  ' ,KS, '  KS/D84  -  ' ,DKS, 

1’  KS/D65  -  ' ,DKSA 

WRITE(12,*)*FG  -  1 ,FG,'  FGP  -  ’ ,  FGP ,  ’  S  -  \S, 

1’  F  -  ’.F,’  X  -  ’.X 
C 

C********************************************************** 

C*  COMPUTE  CS  FOR  PSI  COMPUTATIONS 

C********************************************************** 

c 

D50G-(D65*D35)**0.5 
IF(CS) 125,120,125 
120  CS-D50G*3. 12-0.62 

IF(D50G.LT.0.20)  CS-0.0 
IF(D50G.GT.0. 52)  CS-1.0 

WRITE (12,*) ’ (D65*D35)**0. 5  -’.DSOG,’  CS  -’,CS 
C 

0* *********************************** ********************** 

C*  DETERMINE  PSI  FOR  EACH  BED  MATERIAL  GRAIN  SIZE 
C*  FRACTION 

C********************************************************** 

c 

125  BETA-ALOG 10(10.6) 

BETAX 1 -ALOG 1 0 ( 1 0 . 6  *X) 

RS- (USP**2) / 32 . 2 

DELTA- 11.6*VIS /US 

FI— 1 . 65* (BETA/BET AX 1 ) **2/RS 

DO  135  1-2, ND 

F2«DELTA+CS*GM(I) 

IF(GM(I) .GT.F2)  GO  TO  130 


PSI(I)-F1*F2 
GO  TO  135 

130  PSI(I)-F1*GM(I) 

135  CONTINUE 

WRITE(12,270)  (PSI(I) ,1-2,7) 

WRITE(12,275)  (PSI(I) ,I-8,ND) 

270  FORMAT ('PSI( I) ,1-2,7  -  '.6E10.3) 

275  FORMAT('PSI(I) ,1-8,14  -  '.7E10.3) 

C 

C********************************************************** 

C*  DETERMINE  THE  INTENSITY  OF  TRANSPORT  FOR  EACH  GRAIN 
C*  SIZE 

C********************************************************** 

c 

DO  140  1-2, ND 
PHI(I)-FIG10(PSI(D) 

140  CONTINUE 

WRITE(12,285)  (PHI(I) ,1-2,7) 

WRITE (12, 290)  (PHI (I), 1-8, ND) 

285  FORMAT  ('PHI  (I)  ,1-2,7  -  \6E10.3) 

290  FORMAT  (’PHI  (I),  1-8,  ND  -  \7E10.3) 

C 

C********************************************************** 

C*  COMPUTE  THE  BEDLOAD  DISCHARGE  FOR  EACH  GRAIN  SIZE  RANGE 

C********************************************************** 

c 

DO  145  1-2, ND 

145  IBQB(I)-1200.0*FB(I)*GM(I)**(3.0/2.0)*PHI(I) 
WRITE(12,295)  (IBQB(I) ,1-2,7) 

WRITE(12,300)  (IBQB(I) ,I-8,ND) 

295  FORMAT (' IBQB( I) ,1-2,7  -  '.6E10.3) 

300  FORMAT ('IBQB (I) ,1-8, ND  -  '.7E10.3) 

C 

C********************************************************** 

C*  DETERMINE  THE  CONCENTRATION  IN  THE  BEDLOAD  ZONE 

C********************************************************** 

C 

DO  150  1-2, ND 

CB(I)-IBQB(I) / (11 ,6*2.0*GM(I)*USP) 

WRITE (12 ,^05)  I ,CB(I) 

305  FORMATCCBC  ,12,')  -  '.E10.3) 

C 

C********************************************************** 

C*  COMPARE  CA  TO  CB 

C********************************************************** 

c 

CA(I)»FS(I)*CON*6.24E-5 
IF(CA(I) .LE.CB(I))  GO  TO  150 
CB(I)-1.00i*CA(I) 

IBQB(I)-CB(I)*11.6*2.0*GM(I)*USP 
150  CONTINUE 

WRITE(12,310)  (CA(I) ,1-2,7) 

WRITE (12, 3 15)  (CA(I),I-8,ND) 

310  FORMAT (’CA( I),  1-2, 7  -  \6E10.3) 
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315  FORMAT (’CA(I), 1-8, ND  -  ' .7E10.3) 

WRITE(12,320)  (CB(I) ,1-2 , 7) 

WRITE(12,325)  (CB(I) ,1-8, NO) 

320  FORMAT('CB(I) ,1-2,7  -  ' ,6E10.3) 

325  FORMAT ('CB(I) ,1-8 ,ND  -  '.7E10.3) 

WR1TE(12 ,330)  (IBQB(I) ,1-2 ,7) 

WRITE(12,335)  (IBQB(I) ,1-8 ,ND) 

330  FORMAT (' IBQB (I) ,1-2,7  -  ' .6E10.3) 

335  FORMAT (’ IBQB (I)  ,1-8,14  -  \7E10.3) 

C 

C********************************************************** 

C*  DETERMINE  Z  FOR  EACH  SAMPLED  SUSPENDED  SEDIMENT  GRAIN 
C*  SIZE 

C********************************************************** 

c 

DO  155  1-2, ND 

Fl-0.0 

F2-1.0 

IF(CA(I) .LE.0.0)  GO  TO  155 
IF(CB(I) .LE.O.O)  GO  TO  155 
IF(FS(I).LE. 0.001)  GO  TO  155 
F1-ALOG10(CB(I)/CA(I)) 

F2-2. 0*01(1) 

T-KS/ (30.2*X) 

IF(F2.LT.T)F2-T 
F2-DSV/F2-1.0 
F2-ALOG10(F2/ (DSV/YS-1 .0) ) 

Z(I)-F1/F2 
155  Z1(I)-Z(I) 

WRITE (12, 340)  (Z(I) ,1-2,7) 

WRITE ( 1 2 , 34 5)  (Z(I) ,I-8,ND) 

340  FORMATCZ(I)  ,1-2,7  -  ' ,6E10.3) 

345  FORMAT(’Z(I) ,1-8,14  -  ’.7E10.3) 

C 

C* *************** ********** ******************************** 

C*  COMPUTE  THE  REGRESSED  RELATION  Z-A*W**B 

Q* ************************************************** ******* 

c 

DO  160  1-1 ,ND 

Fl-36.0*VIS**2/(53. 13*GM(I)**3.0) 

160  W(I)-((2.0/3.0+FI)**0.5~F1**0.5)*(53. 13*GM(I))**0.5 
WRITE (12, 350)  (W(I), 1-1,7) 

WRITE(12,355)  (W(I) ,I-8,ND) 

350  FORMAT (’W(I),  1-1,7  -  \7E10.3) 

355  FORMAT('W(I) ,1-8,14  -  ’.7E10.3) 

Z(l)-0.0 

CALL  RGRS(ND,W,Z,FS,AL,B,RSQ,STD) 

C 

Q*********************  **************************  *********** 

C*  RECOMPUTE  Z  FOR  EACH  SIZE  FRACTION 

C********************************************************** 


DO  165  I-l.ND 
Z(I)-AL*W(I)**B 
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165  CONTINUE 

WRITE(12,365)  AL,B, (W(I) ,1-1 ,ND) 

365  FORMAT ( ' AL  -  ',E10.3,/,'B  -  '.E10.3,/, 

l'W(I), 1-1,7  -  ’ ,7E10.3 ,/ , 'W(I) ,1-8 , 14  -  ’.7E10.3) 

WRITE (12, 366) (Z(I),I-1,ND) 

366  FORMAT('Z(I) ,1-1,7  -  '.7E10.3,/, 

rz(l),l-8,14  -  '  ,7E10.3) 

C 

c********************************************************** 

C*  RECOMPUTE  CB  USING  THE  NEW  Z  VALUES  FOR  THE  SAMPLED 
C*  SUSPENDED  SEDIMENT  GRAIN  SIZES 

C********************************************************** 

c 

CA( 1)-FS { 1) *CON*6 . 24E-5 
DO  170  I-l.ND 
IF(FS(I) .IE.0.0)  GO  TO  170 
F1-2.0*GM(I) 

IF(F1.LT.T)F1-T 
F1-DSV/F1-1 .0 
F2-YS/ (DSV-YS) 

CB(I)-CA(I)*(F1*F2)**Z(I) 

170  CONTINUE 
C 

C********************************************************** 

C*  RECOMPUTF,  IBQB  FOR  EACH  SAMPLED  SIZE  FRACTION 

C ********************************************************** 

C 

DO  175  I-l.ND 

175  IBQB(I)-CB(I)*11 . 6*2.0*GM(I)*USP 
WRITE(12,370)  (CB(I) ,1-1 ,ND) 

370  FORMAT(’CBU)  ,1-1,7  -  ' .7E10.3,/, 

1 'CB(I) , 1-8, 14  -  ' , 7E10. 3) 

WRITE(12,371)  (IBQB(I) , I-l.ND) 

371  FORMAT (* IBQB (I) ,1-1 ,7  -  '.7E10.3,/, 

1’ IBQB (I) ,1-8,14  -  T , 7E10. 3) 

C 

C* ******* ************************************************** 

C*  COMPUTE  ISQS  AND  ITQT  FOR  EACH  SIZE  FRACTION 

C********************************************************** 

C 

P-2.303*ALOG10(30. 2*DSV*X/KS) 

WRITE(12,375)  P 

375  FORMAT ( * P  -  ' ,E10.3, // ,3X, ’I’ ,5X, ’ Z(I) ’ ,6X, ' A’ ,9X, 
l'XIl' ,7X,'XI2' ,5X,'ISQS(I)’ ,3X, 'IBQB(I) ’ ,3X, ’ITQT(I) ' ) 
WRITE (*,*)' COMPUTING  CONCENTRATION  PROFILE  INTEGRALS' 
CONV-1 . OE-3 

WRITE (* ,*)  ' CONV  -  ' , CONV 

DO  180  I-l.ND 

ZI-Z(I) 

AI-2.0*GM(I) 

IF(AI.LT.T)AI-T 

A(I)-AI/DSV 

AI-A(I) 

WRITE(* ,*)  ’I  -  ',1 


C  TEST  TO  AVOID  -  REAL  math  overflow  error 
C 

ZMAX-1 .0/ ( .8587E-2-. 1395E-l*ALOG(AI) ) 

XI1-0.0 

XI2-0.0 

IF(ZI.LT.ZMAX)  CALL  I1I2(ZI ,AI ,C0NV,XI1 ,XI2) 

I SQS ( I) -IBQB ( I ) * ( P*XI 1+XI 2 ) 

ITQT(I)-ISQS(I)+IBQB(I) 

180  WRITE (1 2, 380)  I,Z(I) ,A(I) ,XI1 ,XI2,ISQS(I) ,IBQB(I) , 
&ITQT(I) 

380  FORMAT (14 , 7E10 . 3) 

C 

C* *************************************************** ****** 

C*  COMPUTE  THE  TOTAL  LOAD  FOR  THE  SAMPLING  STATION 

C********************************************************** 

C 

TL-0.0 

TG-0.0 

TF-0.0 

TS-0.0 

TFS=0 . 0 

DO  185  I-l.ND 

ITQT(I)-ITQT(I)*43.2*WD 

TL-TL+ITQT(I) 

F1-GM(I)*304.8 

IFCFl.LT. 0.0625)TF-TF+ITQT(I) 

IF (FI . GE . 0 . 0625 . AND . FI . LE . 2 . 0)TS-TS+ITQT (I) 

IF ( F 1 . GT . 2 . 0) TG-TG+ITQT (I ) 

185  CONTINUE 
TFS-TS+TF 

WRITE(12,390)  (ITQT(I) ,1-1 fND) 

390  FORMAT(/,'ITQT(I), 1-1,7  -  ’.7E10.3,/, 

1 ' ITQT(I) ,1*8,14  -  ’.7E10.3) 

WRITE(12,391)TL 

391  FORMATC//TL  -  ’,E10.3,'  TONS/DAY') 

WRITE(13,395)  IST,Q,TL,TF,TS,TG,TFS,X1 ,X2fISTA,D50G,^S 
395  FORMAT(1X,I3,6F10.2,2F6.2,I4,2E10.3) 

900  CONTINUE 
CLOSE (11) 

CL0SE(12) 

CLOSE (13) 

STOP 

END 

C 

C 

FUNCTION  FX(SIN) 

IFCSIN.GE.10.)  THEN 
FX-1 . 

ELSEIF(SIN.LT.O. 16)  THEN 
FX-3.47*SIN 

ELSEIF(SIN.LT.O. 35)  THEN 

FX-0 . 5506*ALOG10(SIN) **2+2 . 42708*ALOG10(SIN)+2 .13114 
ELSEIF(SIN.LT.0.76)  THEN 


FX-1 . 66704+0 . 40926*AL0G10(SIN) -1 . 64269*ALOG10(SIN) **2 
ELSEIF (SIN .LT .0.87)  THEN 

FX-1. 58819-0. 90491*ALOG10(SIN)-7.11837*ALOG10(SIN)**2 
ELSEIF( SIN. LT. 1.20)  THEN 

FX-1 .61092-0. 1472 i*ALOG10 (SIN) -0.80425*ALOG10( SIN) **2 
ELSEIF(SIN.LT. 1.38)  THEN 

FX-1. 56414+1. 02233*ALOG10(SIN)-8.11388*ALOG10(SIN)**2 
ELSEIF(SIN.LT.2.88)  THEN 

FX-1 . 72783-1 . 31612*ALOG10(SIN)+0. 23774*ALOG10(SIN)**2 
ELSEIF ( SIN . LT . 3 . 3 1 )  THEN 

FX-2 . 83698-6 . 13852*AL0G10 (SIN) +5 . 47947*AL0G10(SIN) **2 
ELSEIF(SIN.LT.6.92)  THEN 

FX-1 . 41244-0 . 65949*ALOG10 ( SIN)  -K) . 21 1 18*AL0G10 (SIN) **2 
ELSEIF(SIN.LT. 8. 32)  THEN 

FX-3. 45772-5. 52922*AL0G10(SIN)+3 . 10983*ALOG10(SIN)**2 
ELSE 

FX-5.49124*ALOG10(SIN)-2.87955*AL0G10(SIN) **2-1. 61169 
ENDIF 
RETURN 
END 
C 
C 

SUBROUTINE  FIG5  (X,Y) 

c********************************************************** 

C  THIS  SUBROUTINE  APPROXIMATES  EINSTEINS  FIG  5  BY  A  SERIES 
C  OF  EQNS. 

C ********************************************************** 
IF  (X.LE.1.0)  THEN 
Y-40.0*X**(-1 . 288) 

ELSEIF (X.LE. 2.0)  THEN 
Y-40.0*X** (-0.982) 

ELSEIF (X.LE. 4.0)  THEN 
Y-31 . 1*X** (-0.618) 

ELSEIF (X.LE. 8.0)  THEN 
Y-26.0*X** (-0.486) 

ELSE 

Y-21 . 4*X** (-0. 394) 

ENDIF 

RETURN 

END 

C 

C 

FUNCTION  FIGIO(X) 

C********************************************************** 

C  THIS  FUNCTION  APPROXIMATES  EINSTEINS  FIG  10  BY  A  SERIES 
C  OF  EQNS. 

C********************************************************** 

IF(X.LT.0.4)  THEN 

FIG10-8. 106078*X**(-0. 989002) 

ELSEIF(X.LT.l.O)  THEN 
FIG10-7 . 473444*X** (-1 . 067367) 

ELSEIF (X.LT. 2.0)  THEN 

FIG 10-7 . 512081*X** (-1 . 182832) 

ELSEIF (X.LT. 3.0)  THEN 
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FIG10-9 . 151968*X**(-1 . 472483) 
ELSEIF(X.LT.5.0)  THEN 
FIG10«16.025332*X** (-1.972806) 
ELSE1F(X.LT.7.0)  THEN 
FIG10-41.421873*X**(-2. 578713) 
ELSEIF(X.LT.IO.O)  THEN 
FIG10«l30.2209*X**(-3. 165102) 
ELSEIF(X.LT. 13.0)  THEN 
FIG10-584 . 560594*X** (-3 . 818249) 
ELSEIF(X.LT. 17.0)  THEN 
FIG10«3543.989886*X**(-4. 522935) 
ELSEIF(X.LT. 20.0)  THEN 
FIG10-1690645.945*X** (-6. 703904) 
ELSEIF(X.LT.23.0)  THEN 
FIG10-5188257832.0*X**(-9. 381815) 
ELSE 

FIG10-3 . 134987996E1 5*X** (-13 . 597547) 

ENDIF 

RETURN 

END 


FUNCTION  VISCOS(TEMP) 

DIMENSION  V(13) 

DATA  V/ 2 . 22E-5 , 1 . 92  2E-5 , 1 . 686E-5 , 1 . 49 1 E-5 , 

1 1 . 33E-5 , 1 . 194E-5 , 1 . 08E-5 , 9 . 829E-6 , 8 . 993E-6 , 
28 . 269E-6 , 7 . 636E-6 , 7 . 08E-6 , 6 . 589E-6 / 

TEP- (TEMP-32.0) *5. 0/9.0 
NP0S-TEP/4.0+2 
TO-FLOAT (NPOS) *4 .0-8.0 
IF(TO-TEP)  570,560,560 
560  VISCOS-V(NPOS) 

RETURN 

570  P-(TEP-T0)/4.0 

F1«V(NPOS-1)*(P-P*P)*(P-2.0)/6.0 

F2»V(NPOS) * (P*P-1 . 0) * (P-2 . 0) /2 . 0 

F3-V (NPOS+1 ) * ( 2 . 0*P-P*P) * (P+1 . 0) 12 . 0 

F4-V(NPOS+2) * (P*P*P-P) /6 . 0 

VISCOS-F1+F2+F3+F4 

RETURN 

END 


SUBROUTINE  SHEAR(U,R,RKS,VIS,XF,UST) 

SL-2.0 

SLH-SL 

SLL-0.0 

105  SL* ( SLH+SLL) / 2 . 0 
UST-SQRT(R*SL*32 . 2) 

DELTA- 11.6*VIS/UST 
SIN-RKS /DELTA 
XF-FX(SIN) 

UB-5 . 75*ALOG10 ( 12 . 27*R*XF/RKS) *UST 
CNV-ABS (UB-U) /U 
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IF(CNV.LE. 0.001) GO  TO  110 
IF (DB . GT .0) SLH-SL 
IF (OB . LT . U) SLL-SL 
IF ( SLL . GE . SLH ) SLH-SLH+2 . 0 
GO  TO  105 
110  RETURN 
END 

C 

C 

SUBROUTINE  RGRS (ND ,  X,  Y ,  F ,  A,  B ,  RSQ , STD) 

DIMENSION  X(15) ,Y(15) ,F(15) 

A-1.0 

B-0.0 

RSQ- 1.0 

NP-ND-1 

NP1-NP 

SLNXY  -  0.0 

SLNX  -  0.0 

SLNY  -  0.0 

SLNXS  -  0.0 

SLNYS  -  0.0 

AY  -  0.0 

DO  15  1-2, ND 

IF(X(I) .LE.0.0)  GO  TO  10 

IF(Y(I) .LE.0.0)  GO  TO  10 

IF(F(I).LE. 0.001)  GO  TO  10 

IC-IFIX( 1000 . 0*F(I) ) 

DO  5  IY-l.IC 
AY-AY+Y (I) 

SLNXY  -  SLNXY+ALOG(X(I))*ALOG(Y(I)) 

SLNX  -  SLNX+ALOG (X ( I ) ) 

SLNY  -  SLNY+ALOG (Y (I) ) 

SLNXS  -  SLNXS+ALOG(X(I))**2 
5  SLNYS  -  SLNYS+ALOG (Y(I))**2 
NP-NP+IC-1 
GO  TO  15 
10  NP-NP-1 
NP1-NP1-1 
15  CONTINUE 

IF(NPl-2)  65,40,20 

C 

C********************************************************** 

C*  COMPUTE  A  &  B  WITH  3  OR  MORE  DATA  POINTS 

C********************************************************** 

0 

20  B- SLNXY -SLNX* SLNY /FLO AT (NP) 

B-B/ (SLNXS- (SLNX* SLNX) /FLOAT (NP) ) 

AY-AY /FLOAT (NP) 

SLNYA-SLNY /FLOAT (NP ) 

IF(B. LE.0.0)  GO  TO  22 

A-EXP ( SLNY A-B* SLNX / FLOAT (NP) ) 

RSQ- ( SLNXY- SLNX* SLNY / FLOAT (NP ) ) **2 
RSQ-RSQ/ (SLNXS-SLNX**2/FLOAT(NP) ) 

RSQ-RSQ / ( SLNY  S-SLNY**2 /FLOAT (NP) ) 


DO  35  1-2, ND 
IF(X(I) .LE.0.0)  GO  TO  30 
IF(Y(I).LE.0.0)  GO  TO  30 
IF(F(I).LE. 0.001)  GO  TO  30 
IC-IFIX(1000.0*F(I)) 

YI-A*X(I)**B 
DO  25  IY-1 ,IC 

DSLNYI  -  DSLNYI+(ALOG(Y (I) ) -ALOG(YI) ) **2 
25  DSLNY  -  DSLNY+(Y (I) -AY) **2 
NP-NP+1C-1 
GO  TO  35 
30  NP-NP-1 
35  CONTINUE 

STD- (DSLNYI /FLOAT (NP) ) **0 . 5 

IF(B.GT.O.O)  RETURN 

A-AY 

B-0.0 

RSQ-0.0 

STD- (DSLNY /FLOAT (NP ) ) **0 . 5 
RETURN 
C 

Q* ************************************************* ******** 

C*  COMPUTE  A  &  B  USING  2  DATA  POINTS 

C***********************************«********************** 

c 

40  Yl-0.0 
Y2-0.0 

DO  60  1-2, ND 
IF (Y (I) ) 60,60,45 
45  IF(Y1)50,50,55 
50  Yl-Y(I) 

Xl-X(I) 

GO  TO  60 
55  Y2-Y(I) 

X2-X(I) 

60  CONTINUE 

B-ALOG10(Y1/Y2)/ALOG10(X1/X2) 

A-Yl/X1**B 
IF(B.GT.O.O)  RETURN 
B-0.0 

A-(Yl+Y2)/2.0 

RETURN 


C********************************************************** 

C*  COMPUTE  A  &  B  USING  ONE  DATA  POINT  ASSUMMING 
C*  Y (I)  -  A  *  X(I)  **  0.7 

C********************************************************** 

c 

65  B-0.7 

DO  75  1-2, ND 
IF(Y(I))75,75,70 
70  A-Y (I) /X(I)**0. 7 
75  CONTINUE 


o  o 


112 


RETURN 

END 


SUBROUTINE  IlI2(Z,A,CONV,XIl ,XI2) 

C ********************************************************** 

C*  INTEGRATION  OF  II  AND  12  USING  SIMPSON'S  RULE 

C***************«****************************************** 

C 

WRITE (*,*)  ' Z,A  -  ' ,Z,A 
XI1-0.0 
XI2-0.0 
FI1-0.0 
FI2-0.0 
C 

N-80 

C-.216*A**(Z-1.0)/(1.0-A)**Z 

C 

10  XI1-((1.0-A)/A)**Z 
XI2-XI1*AL0G(A) 

C 

K-l 

DO  60  I-l.N-l 
IF(K-3)  30,20,20 
20  K-l 

30  AN-A+FLOAT ( I) * ( 1 . 0-A) /FLOAT (N) 

FNC- ( ( 1 . O-AN) / AN) **Z 
IF(K-2)  40,50,50 
40  XI1-XI1+4.0*FNC 

XI2-XI2+4 . 0*FNC*ALOG (AN) 

GO  TO  60 

50  XII -XI 1+2.0* FNC 

XI 2 -XI 2+2 . 0*FNC*ALOG(AN) 

60  K-K+l 
C 

XI 1- ( 1 . O-A) / (FLOAT(N) *3 . 0) *XI 1 
XI2- ( 1 . 0-A) / (FLOAT (N) *3 . 0) *XI 2 
CIl-ABS(XIl-FIl) /XII 
CI2-ABS(XI2-FI2)/XI2 

IF(CI1.LE.C0NV.AND.CI2.LE.C0NV)  GO  TO  70 
IF(N.GE.1.0E5)GO  TO  70 
N-N+N 
FI1-XI1 
FI2-XI2 
GO  TO  10 
70  XI1-C*XI1 
XI2-C*XI2 
RETURN 
END 


APPENDIX  C 


FIVEMILE  CREEK  NEAR 
SHOSHONI,  WYOMING 

This  appendix  presents  figures  used  in  this  study  to  estimate  the 
channel  width  and  cross-sectional  area  for  Flvemlle  Creek  near  Shoshoni, 
Wyoming.  The  data  for  these  figures  are  available  in  the  files  of  the 
US  Geological  Survey,  Cheyenne,  Wyoming.  The  data  were  collected 
between  September,  1953  and  October,  1954. 
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APPENDIX  D 


GAGING  DATA  FOR  GOODWIN  CREEK,  MISSISSIPPI 
The  data  in  this  appendix  were  collected  and  provided  by  the 
US  Department  of  Agriculture-Agricultural  Research  Service,  Sedimenta¬ 
tion  Laboratory,  Oxford,  Mississippi.  The  data  were  collected  at  the 
Goodwin  Creek  research  catchment.  Station  2.  The  streamflow  and 
suspended-sediment  data  were  collected  on  August  27,  1986.  The  bed 
material  was  sampled  on  May  9,  1986. 


Table  D.l.  Bed  Material  at  Station  2, 
Goodwin  Creek,  Mississippi 


Size,  mm 

Percent  Finer 

0.0625 

0.24 

0.088 

0.26 

0.125 

0.34 

0.177 

0.74 

0.25 

2.9 

0.354 

6.1 

0.5 

8.0 

0.707 

9.7 

1.0 

12.4 

2.0 

19.6 

4.0 

30.0 

8.0 

46.3 

16.0 

69.6 

38.05 

100.0 
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Table  D.2. 

Streamflov  Data  for  Goodwin  Creek, 
Water  Temperature  72°  F 

Station  2 

Time 

Discharge 

Effective 

Width 

Cross- 

Sectional 

Area 

Hydraulic 

Radius 

Depth  of  Sampled 
Vertical 

hrs 

cfs 

ft 

ft 

ft 

ft 

1935 

440.28 

15.87 

50.46 

1.61 

3.18 

1941 

659.43 

18.46 

69.95 

2.05 

3.79 

1943 

683.71 

18.69 

71.94 

2.09 

3.85 

1944 

725.28 

19.06 

75.30 

2.16 

3.95 

1946 

768.25 

19.43 

78.71 

2.23 

4.05 

1951 

858.40 

20.15 

85.62 

2.37 

4.25 

1954 

905.61 

20.49 

89.14 

2.43 

4.35 

1955 

929.75 

20.66 

90.92 

2.47 

4.40 

1957 

964.15 

20.90 

93.42 

2.51 

4.47 

2001 

989.16 

21.07 

95.22 

2.55 

4.52 

2002 

1019.65 

21.26 

97.39 

2.59 

4.58 

2003 

1014.53 

21.23 

97.03 

2.58 

4.57 

2004 

1014.53 

21.23 

97.03 

2.58 

4.57 

2005 

1014.53 

21.23 

97.03 

2.58 

4.57 

2006 

1019.65 

21.26 

97.39 

2.59 

4.58 

2008 

1014.53 

21.33 

98.12 

2.60 

4.60 

2011 

1035.09 

21.36 

98.48 

2.61 

4.61 

2012 

1040.27 

21.40 

98.85 

2.61 

4.62 

2013 

1045.46 

21.43 

99.21 

2.62 

4.63 

2016 

1035.09 

21.36 

98.48 

2.61 

4.61 

2018 

1035.09 

21.36 

98.48 

2.61 

4.61 

2019 

1040.27 

21.40 

98.85 

2.61 

4.62 

2020 

1035.09 

21.36 

98.48 

2.61 

4.61 

2021 

1035.09 

21.36 

98.48 

2.61 

4.61 

2023 

1029.93 

21.33 

98.12 

2.60 

4.60 

2024 

1029.93 

21.33 

98.12 

2.60 

4.60 

2025 

1019.65 

21.26 

97.39 

2.59 

4.58 

2028 

1014.53 

21.23 

97.03 

2.58 

4.57 

2029 

1009.43 

21.20 

96.67 

2.57 

4.56 

2031 

999.27 

21.13 

95.94 

2.56 

4.54 

2032 

979.11 

21.00 

94.50 

2.53 

4.50 

2033 

979.11 

21.00 

94.50 

2.53 

4.50 

2036 

949.32 

20.80 

92.35 

2.49 

4.44 

2037 

929.75 

20.66 

90.92 

2.47 

4.40 

2040 

881.83 

20.32 

87.38 

2.40 

4.30 

2041 

872.41 

20.25 

86.68 

2.39 

4.28 

2042 

867.72 

20.22 

86.33 

2.38 

4.27 

2043 

849.13 

20.08 

84.93 

2.35 

4.23 

2044 

844.52 

20.04 

84.58 

2.35 

4.22 

2050 

781.41 

19.54 

79.73 

2.25 

4.08 

2056 

704.32 

18.88 

73.62 

2.13 

3.90 

2100 

663.44 

18.49 

70.28 

2.06 

3.80 

(Continued) 
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Time 

hrs 

2105 
2110 
2115 
2110 
2125 
2130 
2135 
2145 
2155 
2205 


Table  D.2  (Concluded) 


Discharge 

cfs 

Effective 

Width 

ft 

Cross- 

Sectional 

Area 

ft 

Hydraulic 

Radius 

ft 

Depth  of  Sampled 
Vertical 
ft 

597.10 

17.82 

64.69 

1.94 

3.63 

534.68 

17.12 

59.22 

1.81 

3.46 

479.47 

16.42 

54.18 

1.70 

3.30 

462.91 

16.19 

52.62 

1.66 

3.25 

385.01 

15.00 

45.00 

1.47 

3.00 

344.90 

14.30 

40.90 

1.40 

2.86 

312.49 

13.70 

37.54 

1.34 

2.74 

250.72 

12.45 

31.00 

1.22 

2.49 

196.52 

11.20 

25.09 

1.10 

2.24 

158.46 

10.20 

20.81 

1.00 

2.04 

Table  D.3.  Point  Sampler  Data  for  Goodwin  Creek,  Station 
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Table  D.4.  Total  Load  Station  Data  for  Goodwin  Creek,  Station 
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